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ABSTRACT 


Results  are  presented  from,  a  group  of  theoretical  studies  of  wave 
propagation  in  horizontally  stratified  media  representative  of  the  ionosphere 
and  the  sub- ionosphere.  The  report  contains:  four  studies  of  the  transient 
three-dimensional  magnetic  fields  produced  by  magnetic  dipoles  in  various 
combinations  of  horizontally  stratified  isotropic  plasma  and  vacuum  regions; 
and  six  studies  of  one- dimensional  hydromagnetic  wave  propagation  in  hori¬ 
zontally  stratified  plasma  and  vacuum  regions.  Also  included  are  biblio¬ 
graphies  concerning  natural  ULF  magnetotelluric  effects,  TjI»F  lightning 
effects,  reflection  and  refraction  of  hydromagnetic  waves,  and  hydromag¬ 
netic  wave  propagation  in  the  ionosphere. 


TABLE  OF  CONTENTS 


Page 


ABSTRACT  iii 

FOREWORD  iv 

LIST  OF  ILLUSTRATIONS  vii 

LIST  OF  TABLES  ix 

LIST  OF  COMMON  SYMBOLS  x 

SECTION  I  INTRODUCTION  1 

SECTION  II  SUMMARY  OF  THREE-DIMENSIONAL  WAVE 

PROPAGATION  STUDIES  3 

2.  1  INTRODUCTION  3 

2. 2  MAGNETIC  DIPOLE  IN  A  UNIFORM  PLASMA  3 

2.  3  HORIZONTAL  DIPOLE  IN  A  SEMI- INFINITE 

PLASMA  8 

2. 4  VERTICAL  DIPOLE  IN  A  SEMI- INFINITE 

PLASMA  9 

2.  5  HORIZONTAL  DIPOLE  ABOVE  A  VACUUM 

CAVITY  1 1 

SECTION  HI  S  UMMARY  OF  ONE- DIMENSIONAL 

HYDROMAGNETIC  WAVE  STUDIES  19 


3.  1 

INTRODUCTION 

19 

3.2 

ESCAPE  FROM  A  VACUUM  CAVITY  INTO 

A  SEMI- INFINITE  PLASMA 

20 

3.3 

ESCAPE  FROM  A  VACUUM  CAVITY 
THROUGH  A  FINITE  PLASMA  SLAB 

29 

3.4 

ESCAPE  FROM  A  VACUUM  CAVITY  INTO 

AN  EXPONENTIAL  PLASMA  REGION 

35 

3.5 

DESCENT  OF  A  HYDROMAGNETIC  WAVE 
INTO  A  VACUUM  CAVITY 

43 

3.6 

REFLECTION  FROM  EXPONENTIAL 
BARRIERS 

47 

3.7 

OTHER  PROBLEMS 

51 

SECTION  IV  APPLICATION  TO  THE  IONOSPHERE 


53 


v 


TABLE  OF  CONTENTS  (cont) 


SECTION  V  CONCLUDING  REMARKS  55 

APPENDIX  A  THREE-  DIMENSIONAL  ELECTROMAGNETIC 

FIELD  PROBLEMS  57 

A.  1  INTRODUCTION  57 

A.  2  BASIC  EQUATIONS  58 

A.  3  MAGNETIC  DIPOLE  IN  A  UNIFORM 

PLASMA  61 

A.  4  HORIZONTAL  DIPOLE  IN  A  SEMI-INFINITE 

PLASMA  65 

A.  5  VERTICAL  DIPOLE  IN  A  SEMI- INFINITE 

PLASMA  69 

A.  6  HORIZONTAL  DIPOLE  ABOVE  A  VACUUM 

CAVITY  73 

APPENDIX  B  ONE- DIMENSIONAL  HYDROMAG NETIC  WAVE 

PROBLEMS  83 

B.  1  INTRODUCTION  83 

B.  2  BASIC  EQUATIONS  84 

B.  3  ESCAPE  OF  A  MAGNETIC  PULSE  FROM  A 

VACUUM  CAVITY  89 

B.  4  DESCENT  OF  A  HYDROMAGNETIC  WAVE 

INTO  A  VACUUM  CAVITY  97 

B.  5  REFLECTION  FROM  AN  UNLIMITED 

EXPONENTIAL  BARRIER  101 

B.  6  REFLECTION  AND  TRANSMISSION  AT  A 

LIMITED  EXPONENTIAL  BARRIER  103 

B.  7  REFLECTION  AND  TRANSMISSION  AT  AN 

ABRUPT  BARRIER  108 

B.  8  ESCAPE  OF  A  MAGNETIC  PULSE  FROM  A 

VACUUM  CAVITY  II  111 

APPENDIX  C  VARIOUS  INTEGRALS  125 

APPENDIX  D  EVALUATION  OF  INTEGRALS  127 

APPENDIX  E  FREQUENCY- TIME  TRANSFORMS  133 

REFERENCES  139 


BIBLIOGRAPHY 


141 


Page 


Figure 

2.  1 

2.2 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 
2.9 
3.  1 

3.2 

3.3 


3.4 


3.5 


3.  6 


3.  7 


3.8 


LIST  OF  ILLUSTRATIONS 


Transient  Magnetic  Field  for  a  Dipole  in  a  Complex 
Plasma 

Magnetic  Field  at  the  Wave  Front 

Transient  Magnetic  Field  for  a  Dipole  in  an  Ohmic 

Plasma 

Two  Equivalent  Problems 

Magnetic  Fields  for  a  Horizontal  Dipole  in  a  Semi- 
Infinite  Plasma 

Magnetic  Field  for  a  Vertical  Dipole  in  a  Semi- 
Infinite  Flasma 

Magnetic  Field  in  a  Cavity  Below  a  Perfect  Plasma 
Magnetic  Field  in  a  Cavity  Below  an  Ohmic  Plasma 
Peak  Magnetic  Field  in  a  Cavity 
Escape  of  a  Magnetic  Field  from  a  Vacuum  Cavity 
Into  a  Semi- Infinite  Plasma 

Escape  of  a  Magnetic  Field  from  a  Vacuum  Cavity  into 
a  Perfect  Plasma;  Spatial  Variation  for  Several  Times 
Escape  of  a  Magnetic  Pulse  from  a  Vacuum  Cavity 
into  an  Ohmic  Plasma;  Time  Histories  of  the  Mag¬ 
netic  Field. 

Escape  of  a  Magnetic  Pulse  from  a  Vacuum  Cavity 
into  an  Ohmic  Plasma;  Spatial  Variation  of  the  Mag¬ 
netic  Field. 

Time  History  and  Polarization  of  Magnetic  Field 
Escaping  from  a  Vacuum  Cavity  into  a  Hall- Effect 
Plasma. 

Time  History  of  Rate- of- Change  of  Magnetic  Field 
Escaping  from  a  Vacuum  Cavity  into  a  Hall- Effect 
Plasma 

Escape  of  a  Magnetic  Field  from  a  Vacuum  Cavity 
Through  a  Finite  Plasma  Slab. 

Decay  of  the  Magnetic  Field  in  a  Cavity  Below  a  Finite 
Ohmic  Slab 


5 

6 

7 

8 

10 

12 

13 

14 
17 

21 

23 

25 

26 

27 

28 

30 

34 


Vll 


LIST  OF  ILLUSTRATIONS  (cont) 

Figure  Page 

3.9  Escape  of  a  Magnetic  Field  Through  an  Exponential 

Plasma  Region  36 

3.  10  Cavity  Magnetic  Response  to  a  Step  Function  Input  39 

3.  1 1  Frequency  Response  Variations  41 

3.  12  Characteristic  Frequencies  of  the  Nighttime 

Ionosphere  42 

3.  13  Descent  cf  a  Hydromagnetic  Wave  into  a  Vacuum 

Cavity  44 

3.  14  Growth  of  the  Magnetic  Field  in  the  Cavity  45 

3.  15  Reflected  Magnetic  Field  46 

3.  16  Unlimited  Exponential  Barrier  Problem  48 

3.  17  Limited  Exponential  Barrier  Problem  48 

3.  18  Reflected  Wave  Shape  49 

3.  19  Energy  Transmission  Through  an  Exponential  Barrier  50 

A.  1  Horizontal  Dipole  Problem  Geometry  62 

A.  2  Horizontal  Dipole  in  a  Semi- Infinite  Plasma  66 

A.  3  Two  Equivalent  Problems  66 

A.  4  Vertical  Dipole  in  a  Semi- Infinite  Plasma  70 

A.  5  Horizontal  Dipole  Above  a  Vacuum  Cavity  74 

B.  1  Escape  of  a  Magnetic  Field  from  a  Vacuum  Cavity  into 

a  Semi- Infinite  Plasma  90 

B.2  Descent  of  a  Hydromagnetic  Wave  into  a  Vacuum  Cavity  98 

B.  3  Unlimited  Exponential  Barrier  Problem  101 

B.4  Comparison  of  Exact  and  Approximate  Expressions  104 

B.  5  Abrupt  Barrier  Problem  108 

B.  6  Multiple  Layer  Problem  112 


viii 


I 

6 

I 


C 

G 

D 

D 


G 

0 


I 

I 


i 

4 


( 


i 

] 


LIST  OF  TABLES 


Table  No. 

I  Values  of  Constants  for  an  Ohmic  Plasma 

II  Values  of  Constants  for  Exponential  Plasma  Problem 
E.  1  Table  of  Transforms 


32 

37 

136 


LIST  OF  COMMON  SYMBOLS 


A(t) 

B* 


B 

z 


B 


P* 


c 


C 


erf 

erfc 

jir 

E  ,E  ,E 
x  y  z 

f 

F 

h 

h(t) 

ha 

H 

H(iu) 

I 

J,  J' 

J 

J  * 

k 

K 


m 


Response  of  a  linear  system  to  a  unit  step  function 
Time-dependent  component  of  magnetic  induction  vector 
Constant  vector  component  of  a  magnetic  field 
Components  of  B  in  rectangular  cartesian  coordinates 
Components  of  B  in  spherical  polar  coordinates 
Speed  of  light  in  vacuum 
Dipole  moment  times  P-  /4^ 

Cosine  integral 
Error  function 

Complementary  error  function 
Electric  Field  vector 
Rectangular  cartesian  components  of  E 
Frequency  in  cycles  per  second 
Vector  potential 

Height  above  sub-ionosphere  or  vacuum  cavity 
Response  of  a  linear  system  to  a  unit  impulse 
Height  of  plasma  slab  above  a  vacuum  cavity 
Height  of  sub- ionosphere  or  vacuum  cavity 
Frequency  response  function  of  a  linear  system 
Modified  Bessel  function  of  the  first  kind 
Bessel  function  of  the  first  kind  and  its  first  derivative 
Current  density  vector 

Asymptotic  expression  for  J  (See  Eqs.  (D.17)  and  (D.18)) 
Wave  number 

Modified  Bessel  function  of  the  second  kind 
Dipole  moment 


x 


I  1 

ri  » 


sgn  (x) 


Scale  distance  associated  with  Aifven  speed  variation  in 
the  upper  ionosphere 

Radial  distance 

Reflection  coefficient 

Designates  real  part  of  x 

Transformed  \ariable  of  the  Laplace  transform 
Sign  of  x  1  for  x  >  0;  -  1  for  x  <  0) 

Step  function  (S  (t)  =  0  for  t  <  0  and  S  (t)  =  1  for  t  >  0) 

Sine  integral 
Time 


Period  or  characteristic  time 


Plasma  velocity 
Wave  propagation  speed 
Aifven  speed 

'Wave  front  speed  in  a  hydromagnetic  plasma 
Function  defined  by  Equation  (B.  2 .22) 
Function  defined  by  Equation  (B.  8.35) 
Horizontal  distance 

Horizontal  distance  perpendicular  to  X-axis 


Vertical  distance 


Damping  parameter  (Eq.  (2.1)) 

Propagation  constant  (y  =  ik) 

Delta  or  impulse  function  of  t  (see  Appendix  E) 
n’th  derivative  of  the  delta  function  (see  Appendix  E) 


xi 


0 


P 

v 

P 

or 

Or 

O 

7 

<P 

$ 

$ 

* 

(0 


Polar  angle 

Permeability  in  vacuum 
Collision  or  damping  frequency 
Horizontal  distance 
Conductivity 

Real  part  of  conductivity 
Non-dimensional  time 

Acute  angle  between  geomagnetic  field  vector  and  the 
vertical 

Flux 

Divergence  of  vector  potential  in  Appendix  A 
Azimuthal  angle 
Frequency,  rad /sec 


SUBSCRIPTS 

P 

v 


y 

z 


Plasma 
Vacuum 
X-  component 
Y- component 
Z- component 


SECTION  I 


INTRODUCTION 

The  subject  of  ultra- low- frequency  (ULF)  magnetohydrodynamic  wave 
propagation  has  recently  become  of  importance  because  of  its  relationships 
to  solar  induced  geomagnetic  disturbances  and  to  near  1  cps  magnetotelluric 
effects  of  high  altitude  nuclear  detonations.  However,  the  present  under¬ 
standing  of  this  subject  is  quite  limited  because  of  the  considerable  mathe¬ 
matical  difficulties  involved  in  solving  any  problem  of  hydromagnetic  wave 
propagation  through  complex  plasma  media,  except  for  the  simplest  cases 
where  the  plasma  is  completely  uniform  or  lias  slowly  varying  properties  (so 
that  a  WKB  approach  may  be  used).  In  order  to  improve  present  understand¬ 
ing  of  hydromagnetic  wave  propagation  in  the  ionosphere,  two  groups  of  related 
theoretical  wave  propagation  studies  were  made  under  USAECOM  sponsorship, 
the  results  of  which  are  summarized  in  Sections  II  and  in  of  this  report.  The 
first  group  of  studies,  summarized  in  Section  II  and  described  in  detail  in 
Appendix  A,  considers  the  three-dimensional  electromagnetic  fields  produced 
by  magnetic  dipoles  in  various  combinations  of  horizontally  stratified  isotropic 
plasma  regions.  The  second  group  of  studies,  summarized  in  Section  III  and 
described  in  detail  in  Appendix  B  considers  various  one- dimensional  problems 
of  hydromagnetic  wave  propagation.  Section  IV  discusses  the  application  of 
these  study  results  to  ionospheric  wave  propagation  problems. 

A  bibliography  is  included  covering  four  groups  of  topics  associated 
with  the  subject  of  ULF  wave  propagation  in  horizontally  stratified  media 
representative  of  the  ionosphere. 


SECTION  II 


SUMMARY  OF  THREE-DIMENSIONAL  WAVE  PROPAGATION 

STUDIES 

2.  1  INTRODUCTION 

Appendix  A  of  this  report  presents  a  theoretical  analysis  of  four  three- 
dimensional  problems  concerned  with  the  steady-state  and  transient  magnetic 
fields  produced  by  magnetic  dipoles  in  various  types  of  horizontally  stratified 
plasma  and  un-ionized  media.  The  principal  results  of  these  studies  are  sum¬ 
marized  below,  details  being  given  in  Appendix  A. 

It  is  assumed  throughout  this  section  that  Maxwell's  equations  are 
applicable  with  an  isotropic  conductivity  law  and  that  the  displacement  current 
may  be  neglected  in  vacuum  regions.  This  assumption  of  isotropic  conduc¬ 
tivity  is,  of  course,  not  generally  exactly  applicable  to  hydromagnetic  wave 
propagation  studies,  which  are  our  principal  concern  in  this  report.  However, 
to  a  considerable  degree,  both  qualitatively  and  quantitatively,  a  study  of  these 
problems  offers  considerable  insight  into  the  corresponding  hydromagnetic 
problems. 

In  this  section  the  term  perfect  plasma  is  used  to  refer  to  a  plasma 
with  purely  imaginary  conductivity;  the  term  complex  plasma  for  a  complex 
conductivity;  and  the  term  ohmic  plasma  for  a  real  conductivity  (see  Sec. 

A.  2). 

2.  2  MAGNETIC  DIPOLE  IN  A  UNIFORM  PLASMA 

Section  A.  3  of  Appendix  A  presents  a  derivation  of  the  steady-state 
and  step-function-response  electromagnetic  fields  produced  by  a  magnetic 
dipole  in  an  infinite  uniform  plasma.  These  results  are  exactly  applicable  to 
the  analogous  case  of  a  hydrom'gnetic  plasma,  provided  that  the  dipole  axis 
is  parallel  to  the  ambient  magnetic  field  for  the  hydromagnetic  case. 


Considering  only  the  magnetic  fields  produced  by  a  suddenly  created 
magnetic  dipole,  the  transient  magnetic  fields  ara  shown  in  Figures  2.  1 
through  2.3,  where  the  following  symbols  are  used 

c  is  a  damping  parameter 

r  is  distance  from  dipole  to  observer 

v  is  a  characteristic  collision  damping  frequency 

cr^  is  the  real  component  of  conductivity 

V  is  the  phase  speed 

P  is  the  permeability  in  vacuum 

a=4-vr/V  =  j^qVi  (2.  1) 

ana  where  Br  and  B  ^  are  radial  and  tangential  components  of  the  transient 
magnetic  field  in  polar  coordinates  with  respect  to  the  dipole  axis  in  a  plane 
through  the  dipole  axis.  The  azimuthal  component  of  the  transient  magnetic 
field  is  zero. 

For  a  perfect  plasma,  the  transient  magnetic  field  takes  the  form  of 
a  wave  front  moving  into  the  plasma  at  the  phase  speed  V,  with  the  magnetic 
field  ahead  of  the  front  being  zero  (for  Vt/r  <  1)  and  behind  the  front  being 
equal  to  that  for  the  static  dipole  condition  (see  curve  for  a  =  0  in  Fig.  2.  1). 
Both  of  the  magnetic  field  components  B^  and  B^  contain  infinite  impulses 
at  the  wave  front  (see  Eqs.  (A.  3.  13)  and  (A.  3.  14)). 

For  the  case  of  a  complex  plasma  ( c  >  0)  with  small  to  moderate 
value  of  damping  (a  £  1) ,  a  similar  behavior  occurs  (see  curves  for  a  -  1  in 
Fig.  2.  1),  except  that  the  initial  peak  field  strength  (just  behind  the  wave 
front)  for  the  radial  component  (Br)  is  somewhat  smaller  than  the  static  value 
and  is  larger  for  the  tangential  component.  Both  components  vary  monot- 
onically  toward  the  static  value. 

For  greater  values  of  damping  ratios,  the  magnetic  field  at  the  wave 
front  varies  as  indicated  in  Figure  2.  2,  according  to  Equations  (A.  3„  7)  and 
(A.  3.  8).  The  radial  component  decreases  monotonically  with  increased 
damping;  the  tangential  component  increases  with  increased  damping  up  to  a 
damping  ratio  of  about  a  =  3,  and  decreases  thereafter. 

Time  histories  of  the  magnetic  field  for  a  large  value  of  damping  are 
shown  in  Figure  2.  1  for  a  =  5.  It  may  be  noted  that  the  tangential  magnetic 
tield  increases  above  the  initial  values  at  the  wave  front  to  a  maximum  before 
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approaching  the  static  condition;  the  radial  component  increases  monot- 
onically  from  the  initial  value  toward  the  static  value. 

Considering  finally  the  case  of  an  ohmic  plasma,  Figure  2.  3  presents 
time  histories  of  the  corresponding  magnetic  field  components.  These  curves 
correspond  to  the  complex  case  where  a  -*■  oo  and  are  similar  to  the  above 
discussed  curves  for  a  -  5,  except  that  the  initial  value  is  zero  for  the  ohmic 
case,  there  being  no  discinct  finite -amplitude  wave  front  in  this  case. 

2.  3  HORIZONTAL  MAGNETIC  DiPOLE  IN  A  SEMI-INFINITE  PLASMA 

Section  A.  1  of  Appendix  A  considers  the  transient  magnetic  fields 
produced  by  a  horizontal  magnetic  dipole  located  at  a  height  h  above  a  hor¬ 
izontal  interface  separating  a  plasma  above  from  a  vacuum  below. 

The  following  reciprocity  law,  indicated  in  Figure  2.4,  is  found  to  be 

o's  Designate  dipole  locations 
x's  Designate  observation  points 


Vacuum 

(a.)  Reference  Problem  (b.)  Reciprocal  Problem 


Fig. 2. 4  Two  Equivalent  Problems 
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valid  in  the  plasma.  The  magnetic  field  at  tbv.  plasma-vacuum  interface 
created  by  a  dipole  at  the  height  h  (Fig.  2.4a)  is  the  same  as  the  magnetic 
field  created  at  the  height  h  by  a  dipole  at  the  interface  {Fig.  2,4b). 

For  the  special  case  where  both  dipole  and  observation  points  are 

located  in  the  plane  of  the  interface,  the  transient  magnetic  fields  produced 

by  a  suddenly  created  (step-function)  magnetic  dipole  vary  as  indicated  in 

Figure  2.  5.  This  figure  presents  time  histories  of  the  horizontal  magnetic 

field  components  for  a  perfect  plasma  and  an  ohmic  plasma,  where  B  and 

P 

Bg  are  the  radial  and  tangential  components  of  the  horizontal  magnetic 
field. 

For  both  perfect  and  ohmic  plasmas,  the  magnetic  field  of  the  dipole 
immediately  fills  the  vacuum  half-space  {since  displacement  current  is 
neglected),  but  requires  a  finite  time  to  penetrate  into  the  plasma.  Hence, 
the  plasma  region  initially  acts  as  a  perfect  reflector,  creating  twice  the 
static  magnetic  field  at  time  zero  in  the  vacuum  region  (and  at  the  interface). 
With  increasing  time,  the  doubled  magnetic  field  at  the  interface  then  begins 
to  decrease  as  the  excess  dipole  field  penetrates  into  the  plasma.  For  a 
perfect  plasma,  the  magnetic  field  components  eventually  reverse  their 
polarity,  reaching  a  large  negative  peak  at  the  time  t  =  r/V  (Fig.  2.  fa), 
which  is  the  time  required  for  the  dipole  disturbance  to  travel  in  the  plasma 
from  the  dipole  to  an  observer  point.  The  magnetic  field  components  then 
abruptly  reverse  their  direction  and  reach  their  static  values.  For  the  case 
of  an  ohmic  plasma,  the  transition  from  the  initial  doubled  field  to  the  final 
static  state  is  almost  a  simple  monctonic  decay  (Fig.  2.  5b),  with  a  charac¬ 
teristic  time  constant  on  the  order  of  0.  Iuct  r2. 

^  o 

2.  4  VERTICAL  MAGNETIC  DIPOLE  IN  A  SEMI-i  v'FINITE  PLASMA 

Section  A.  5  of  Appendix  A  considers  the  transient  magnetic  fields 
produced  by  a  vertical  magnetic  dipole  located  above  a  horizontal  interface 
separating  a  plasma  above  from  a  vacuum  below.  These  results  are  exactly 
applicable  to  the  analogous  case  of  a  hydromagnetic  plasma,  provided  that 
the  ambient  magnetic  field  vector  for  the  hydromagnetic  plasma  is  a  vertical 
vector. 
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Fig.  Z.  5  Magnetic  Fields  for  a  Horizontal  Dipole  in  a  Semi-infinite  Plasma 


For  the  special  case  of  a  perfect  plasma  where  both  dipole  and  ob¬ 
servation  points  are  located  in  the  plane  of  the  interface,  the  transient  mag¬ 
netic  fields  produced  by  a  suddenly  created  (step -function)  magnetic  dipole 
vary  as  indicated  in  Figure  2.  6,  where  Bp  and  Bz  are  the  radial  and  vertical 
magnetic  field  components,  the  other  component  being  zero.  The  initial  mag¬ 
netic  field  components  are  zero,  since  the  plasma  region  initially  acts  as  a 
perfect  reflector  (see  Sect.  2.  3),  thereby  creating  an  image  dipole  canceling 
the  applied  dipole.  As  time  increases,  the  magnetic  field  components  increase 
to  large  multiples  of  the  static  dipole  field,  reaching  infinity  at  the  time 
t  =  p  /V,  which  is  the  time  required  for  the  dipole-produced  wave  front  to 
travel  in  the  plasma  from  dipole  to  observer.  For  larger  times,  the  vertical 
component  of  the  magnetic  field  immediately  reaches  the  static  dipole  value 
and  the  horizontal  component  decreases  rapidly  toward  zero. 

2.  5  HORIZONTAL  DIPOLE  ABOVE  A  VACUUM  CAVITY 

Section  A.  6  of  Appendix  A  considers  the  problem  where  a  horizontal 
magnetic  dipole  is  placed  in  a  plasma  at  an  attitude  h  above  a  horizontal 
vacuum  cavity  of  depth  H,  the  lower  edge  of  the  cavity  being  the  surface  of  a 
perfect  conductor.  For  the  case  of  a  suddenly  applied  (step-function)  dipole 
in  the  plasma,  the  transient  magnetic  field  in  the  cavity  at  the  surface  of  the 
perfect  conductor  behaves  as  follows  for  perfect,  dissipative,  and  ohmic 
plasma. 

Time  histories  of  the  transient,  magnetic  fields  at  the  surface  of  the 
perfect  conductor,  according  to  Equation  (A.  6.  29),  are  shown  in  Figures  2.  7 
and  2.  8  for  perfect  and  chmic  plasma,  respectively,  for  the  parameter  ratios: 

h/H  =  1  ;  H/p  =0.1  (A.  6.41) 

where 

p  is  the  horizontal  distance  between  dipole  and 
observation  point 

Considering  first  the  perfect  plasma,  it  may  be  seen  from  Figure  2.  7 
that  the  transient  magnetic  field  is  zero  until  the  time  t  =  h/  V,  which  is  the 
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time  required  for  the  wave  front  to  travel  vertically  downward  from  the  dipole 
altitude  to  the  upper  edge  of  the  vacxuim  cavity.  At  this  time,  the  magnetic 
field  suddenly  jumps  to  the  large  multiple  p  / H  of  the  static  value  of  the  dipole 
field  for  the  0-  component  of  the  magnetic  field  and  to  the  multiple  j-  p/H 
for  the  radial  component.  The  magnetic  field  then  decays  exponentially  toward 
zero  as  the  magnetic  field  in  the  cavity  decays  into  the  plasma  (see  Section 
B.  3  of  Appendix  B).  The  characteristic  decay  time  here  is  H/V.  For  large 
values  of  time,  near  t  =  p/V,  the  magnetic  field  in  the  cavity  rises  again  and, 
after  a  short  transient  period-  reaches  the  static  dipole  field  configuration. 

This  second  disturbance  corresponds  to  direct  travel  of  the  electromagnetic 
wave  through  the  plasma  to  the  observation  point  p,  as  contrasted  to  the  early 
time  peak  at  Vt  =  h  which  corresponds  to  a  more  rapid  travel  via  the  vacuum 
cavity. 

Next,  consider  the  case  of  a  cavity  below  an  ohmic  plasma  (Fig.  2.  8). 

If  the  dipole  is  located  at  the  upper  edge  of  the  cavity  or  in  the  cavity  (h/H  =  0), 
the  initial  magnetic  field  is  the  same  as  for  a  perfect  plasma,  and  subsequently 
decays  slowly  toward  the  static  field  configuration  (Fig.  2.  8b),  There  is  no 
conspicuous  second  peak  or  rapid  transition  here  of  the  type  noted  above  for 
a  perfect  plasma,  because  of  the  large  dissipation  of  transients  experienced 
by  waves  traveling  any  large  distance  inside  an  ohmic  plasma. 

For  dipoles  located  above  the  cavity  in  an  ohmic  plasma,  the  peak 
magnetic  field  strength  decreases  with  increasing  altitude  (see  Fig.  2.8  and 
solid  line  in  Fig.  2.  9)  and  occurs  at  progressively  later  times  as  the  altitude 
h  is  increased. 

For  dipoles  located  in  a  dissipative  plasma  above  a  vacuum  cavity,  the  mag¬ 
netic  field  behavior  will  be  intermediate  between  that  for  a  perfect  and  for  an 
ohmic  plasma.  More  specifically,  there  will  be  an  initial  delay  until  the  wave 
front  travels  from  the  dipole  to  the  upper  surface  of  the  cavity.  The  magnetic 
field  in  the  cavity  will  then  jump  suddenly  to  the  value  given  by  Equation 
(A.  6.  35),  which  is  plotted  as  dashed  lines  in  Figure  2.  9.  If  the  dissipation 
is  large,  the  magnetic  field  will  subsequently  increase  to  a  maximum  and  then 
decay  toward  the  static  value;  otherwise,  for  small  dissipation,  the  trends  of 
Figure  2.  7  for  a  perfect  plasma  will  be  followed.  The  peak  value  of  the  mag¬ 
netic  field  transient  may  be  calculated  from  Equation  (A.  6.  35)  for  small  dis¬ 
sipation  (dashed  lines  in  Fig.  2.  9)  and  may  be  estimated  from  the  solution  for 
an  ohmic  plasma  for  large  values  of  dissipation  (solid  line  in  Fig.  2.9)r 
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SECTION  in 


SUMMARY  OF  ONE- DIMENSIONAL  HYDROMAGNETIC 
WAVE  PROBLEMS 

3.  1  INTRODUCTION 

This  Section  presents  a  summary  of  a  group  of  theoretical  studies, 
presented  in  Appendix  B,  which  are  concerned  with  one  dimensional  hydro- 
magnetic  wave  propagation  in  stratified  plasma  and  vacuum  regions.  These 
studies  are  considered  relevant  to  an  understanding  of  near  1  cps  fluctuations  in 
the  earth's  magnetic  field  in  and  below  the  ionosphere,  which  may  be  initially 
excited  by  hydromagnetic  or  electromagnetic  waves  created  by  distant  high 
altitude  nuclear  detonations.  They  are  also  relevant  to  an  understanding  of 
the  ultra- low- frequency  worldwide  effects  of  lightning  discharges. 

Five  types  of  plasma  and  vacuum  regions  are  considered  in  this  sec¬ 
tion,  which  are  designated  as  perfect  plasma,  complex  plasma,  ohmic  plasma. 
Hall- effect  plasma,  and  vacuum.  Precise  definitions  of  the  properties  of 
these  regions  are  given  by  the  equations  of  Section  B.  2  and  arc  briefly  out¬ 
lined  below. 

A  perfect  plasma  is  defined  herein  as  a  conventional  cold  collision- 
free  hydromagnetic  plasma  without  Hall  effects  present  (Eqs.  (B.  2.  5)  and 
(B.2.6)). 

A  complex  plasma  is  defined  as  a  cold  hydromagnetic  plasma  with 
collision  damping  present,  but  without  Hall  effects  (Eqs.  (B.  2.7)  and  (B.2.  8)). 
This  type  of  plasma  simulates  the  F- layer  of  the  ionosphere. 

An  ohmic  plasma  is  defined  as  a  plasma  whose  ratio  of  current  density 
to  electric  field  intensity  is  given  by  the  scalar  Ohm's  law  expression  for  an 
simple  conductor  (J  =  <r  E)  with  a  real  scalar  conductivity  constant. 

A  Hall- effect  plasma  is  defined  as  a  hydromagnetic  plasma  charac¬ 
terized  by  finite  real  Pedersen  and  Hall  conductivity  constants  (Eq.  (B.2.  9)). 
and  infinite  parallel  conductivity.  This  type  of  plasma  simulates  the  E- 
layer  of  the  ionosphere. 


A  vacuum  is  defined  as  a  non-conducting  region  wherein  Maxwell's 
displacement  current  may  be  neglected.  For  frequencies  ne^.r  or  below 
1  cps,  this  type  of  region  simulates  the  behavior  of  the  atmosphere  and 
ionosphere  up  to  about  90  km  altitude  in  the  daytime,  and  uo  to  about  290  km 
at  night. 

Throughout  this  section  it  is  assumed  that  all  plasma  regions  are 
permeated  by  a  uniform  ambient  geomagnetic  field,  whose  field  vector  Bq 
lies  in  the  YZ-plane  of  a  rectangular  cartesian  coordinate  system  having  a 
horizontal  XY- plane 


3.  2  ESCAPE  FROM  A  VACUUM  CAVITY  INTO  A  SEMI- INFINITE  PLASMA 


3.2.  1 


Function  Response 


Section  B.  3  of  Appendix  B  considers  the  escape  of  a  magnedc  pulse 

from  a  vacuum  cavity  of  height  H,  which  is  bounded  below  by  a  perfect 

conductor  and  above  by  a  semi- infinite  hydromagnetic  plasma  containing  an 

ambient  constant  magnetic  field  vector  B^  (Fig.  3„  1).  It  is  assumed  that  at 

o 

time  zero  a  uniform  small  horizontal  magnetic  field  of  strength  Bvq  and  total 
flux  $  =  is  suddenly  applied  to  the  cavity  and  is  then  permitted  to 

relax  into  the  plasma. 

For  the  case  of  a  perfect  hydromagnetic  plasma  (see  Section  B.  2  of 
Appendix  B),  the  initial  magnetic  field  in  the  cavity  decays  exponentially  with 
increasing  time  t  (see  Eq.  (B.  3.  15))  as 


B(t)  =  B  exp  (-  V  t/H) 
vo  r  m 


(3.  1) 


where 

B  is  the  transient  magnetic  field  strength  in  the  cavity 
H  is  the  cavity  depth 

and  where  Vm  is  the  phase  velocity  in  the  hydromagnetic  plasma,  which  is 
(from  Eq.  (B.2,20)) 


(3.2) 


V  for  modified  Alfven  waves 
a 


V  cos  0  for  Alfven  waves 
a 


where 


is  the  Alfven  speed  (V  =  p) 

is  the  ambient  magnetic  field  strength 

is  the  acute  angle  between  the  ambient  mag¬ 
netic  field  vector  B  and  the  vertical. 

c 


and  where  the  component  of  the  magnetic  field  in  the  plane  farmed  by  the 
ambient  magnetic  field  vector  and  the  vertical  behaves  as  a  modified  Alfven 
wave  and  the  component  perpendicular  to  this  plane  vehaves  as  an  Alfven  wave 
wave  (see  Section  B.2).  The  corresponding  transient  magnetic  field  in  the 
plasma  consists  of  a  shock  front  of  amplitude  Bvq  with  a  spatially  exponen¬ 
tially  decaying  wake  moving  into  the  plasma  from  the  cavity  wall,  traveling 
at  the  phase  speed  V  ,  as  is  shown  in  Figure  3.2,  according  to  Equation 
(B.  3. 15). 

For  a  complex  plasma,  more  complicated  equations  apply  (see  Eq. 

(B.  3.11)),  but  the  behavior  will  be  qualitatively  the  same  for  small  values 
of  the  damping  ratio  vH/V^,  where  v  is  a  characteristic  ion  collision 
frequency.  For  large  damping  ratios  (vH/V  -^1),  the  behavior  approaches 
that  for  a  simple  ohmic  plasma  in  that  the  shock  front  vanishes,  the  magnetic 
field  in  the  cavity  B^  decays  as  (from  Eq.  (B.  3.  21)) 


where 


B  =  B  w  (i /t/T) 
v  vo 


T  =  H-  O'  H2 

w(z)  =  exp  (-  z3)  erfc  (-  iz) 
<r  is  the  conductivity 


(3.3) 


(3.4) 

(3.5) 


Escape  of  a  Magnetic  Field  from  a  Vacuum  Cavity  into  a  Perfect  Plasma; 
Spatial  Distribution  for  Several  Times 


riCT-rrrru 


(see  upper  cum  in  Fig.  3.  3),  and  the  wave  in  the  plasma  takes  the  form 
indicated  in  Figures  3.  3  and  3.4  (See  Eq.  (B.  3.19)). 

It  may  he  noted  here  that  the  characteristic  time  associated  with  the 
decay  o£  the  magnetic  field  in  the  cavity  below  the  perfect  plasma  is  H/V  , 
which  is  of  the  order  of  one  second  for  the  nighttime  sub- ionosphere  cavity 
(e.g.,  for  H-SJOG  km,  V  *'-300  km  per  sec).  The  characteristic  decay  time 
for  an  ohmic  plasma  is  )ur  H2  ,  which  is  of  the  order  of  one  second  for  the  day¬ 
time  sub- ionosphere  cavity  (e.g.,  using  a  Pederson  conductivity  for  the  E- 
layer  of  about  2  X  lO-*  MKS  units  with  an  ionosphere  height  of  80  km). 

For  escape  into  an  Hall  effect  plasma,  the  magnetic  field  in  the  cavity 
obeys  a  relationship  of  the  form  (see  Eq.  (B.  3.  34) 


where 


B^t)  =  Cn  w+  +  Cis  w_ 


Bvy(t)  =  Czi  w-r  *  Caa  w- 


wA  =  w  (i/t/Tj. ) 


(3.  6) 


and  where  the  T^’  s  are  complex  constants  (Eq.  (B.  3.  32))  and  where  the  C’  s 

are  each  linear  functions  of  both  of  the  initial  components  B  and  B  of 

vxo  vyo 

the  magnetic  field  in  the  cavity  (see  Eq.  (B.  3.  34)).  Figures  3.  5  and  3.6 
illustrate  this  magnetic  field  variation  and  its  rate  of  change  for  parameters 
appropriate  to  the  daytime  E-layer  of  the  ionosphere  for  the  case  where  the 
ambient  magnetic  field  is  vertical  (high  latitude  condition)  and  where  the 
initial  magnetic  field  is  polarized  in  the  Y- direction  (B^  =  0).  The  magnetic 
field  in  the  Y- direction  decays  in  about  the  same  maimer  as  for  an  ohmic 
plasma  (Fig.  3,3).  A  substantial  X- component  is  seen  to  come  into  existence 
in  less  than  a  second  and  subsequently  decays  at  a  rate  similar  to  that  for 
the  Y-  component. 


3.  2.3  Frequency  Response 

i-Iext,  consider  the  cavity  frequency  response  to  a  steady  state 
sinusoidal  input  of  magnetic  flux.  For  the  cases  of  a  plasma  without  Hall 
effects,  this  is  given  by  Equation  (B.  8.17)  as 
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ime  History  and  Polarization  of  Magnetic  Field  Escaping  from 
Vacuum  Into  a  Hall- Effect  Plasma 


.xtie  History  of  Rate  of  Change  of  I 
Vacuum  Cavity  Into  a  Hall- Effect 


H(iu)  =  ikH/(l  +  ikH) 


where 

H(iu)  is  the  frequency  response  function 
k  is  the  wave  number  in  the  plasma 


and  from  Equations  (B.  2.  28)  through  (B.  2.  30) 


r 


(complex  plasma) 
(perfect  plasma) 

(ohmic  plasma) 


3.  3  ESCAPE  FROM  A  VACUUM  CAVITY  THROUGH  A 
FINITE  PLASMA  SLAB 


Section  B.  8.  3  of  Appendix  B  considers  the  problem  of  the  escape  of 
a  magnetic  pulse  from  a  vacuum  cavity  of  i.eight  H,  which  is  bounded  below 
by  a  perfect  conductor  (where  E  =  0)  and  is  bounded  above  by  a  finite  height 
dense  plasma  slab  of  height  hs,  and  which  is  in  turn  bounded  above  by  a 
rarefied  plasma  region  (Fig.  3,  7).  Also  considered  is  the  response  of  the 
cavity  to  a  sinusoidal  application  of  magnetic  flux  to  the  cavity. 

3.  3.1  Step  Function  Response 

First  consider  the  case  of  the  step- function  response  of  the  cavity. 

It  is  assumed  that  at  time  zero  a  uniform  small  horizontal  magnetic  field  of 
unit  strength  (Bvq  =  1)  and  total  flux  $  =  BvqH  is  suddenly  applied  to  the 
cavity  and  is  then  permitted  to  relax  into  the  plasma. 

Complex  plasma 

For  the  case  of  a  complex  plasma,  the  magnetic  field  in  the  cavity 


varies  as 


Fig.  3.  7  Escape  of  a  Magnetic  Field  from  a  Vacuum  Cavity 
Through  a  Finite  Plasma  Slab 


(B.8.  30a) 


B^(t)  =  1  -  2SCn(l  -  exp  (-  5vt)[cos  +  sfv.lo^)  sinw^t] ) 

=  ZSC^exp  (-  £vt)[cos  o^t  +  ^(v/u)sin  w^t]  (B.  8.  30b) 

where 

v  is  a  collision  damping  frequency 

and  where  the  C  's  and  w  's  are  numerical  constants  defined  in  Section  B.8.  3 
n  n 

of  Appendix  B.  (Table  I  gives  the  first  few  values  of  2Cn  for  a  =  H/hjj  =  3  and  2.) 
The  response  behavior  consists  of  a  superimposition  of  damped  oscillations  at 
the  characteristic  frequencies  un. 

Perfect  plasma 


For  a  perfect  plasma,  the  magnetic  field  in  the  cavity  varies  accord¬ 
ing  to  Equation  (B.  8.  30)  above  with  v  =  0,  or 

B  (t)  =  2SC  cosm  t  (3.8) 

v  n  n 

where 


«  "  £  „v  /ha  , 

n  n  m 


(3.9) 


where  the  !n's  are  the  positive  roots  of  the  equation 


at,  tan  £  -  1 

n  n 


(3.10) 


(see  Table  I  for  the  first  few  values  of  |  for  a  =  3  and  2),  where 

v  3n 


a  =  H/h2 


(3.11) 


h2  is  the  height  of  the  finite  plasma  slab 
V  is  the  phase  speed  in  the  plasma 

The  response  behavior  is  seen  to  consist  of  a  superimposition  of 
continuous  oscillations  at  the  various  characteristic  frequencies  at  the 
amplitudes  2Cn  (see  Eq.  (3.8);  also  Table  I  for  the  first  four  values  of  2Cn 
for  a  =  3  and  2), 
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H/ha 

n 

| 

2C 

T  (sec) 

n 

n 

ni 

1 

.9$ 

.  635 

2.900 

2 

3.54 

.  184 

.226 

2/3 

3 

6.51 

.065 

.067 

4 

9.  58 

.031 

.031 

1 

.65 

.851 

0.  740 

2 

2 

3.29 

.086 

.029 

3 

6.36 

.  024 

.  008 

*  Listed  values  of  £  and  2C  also  apply  for 

n  n 

complex  and  perfect  plasmas 
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Tnx  =  li<rHa/(aS|na)  (3‘13) 

and  where  the  first  few  values  of  C  and  £  are  given  in  Table  I  for 

n  n 

a  =  H/hg  =  5  and  2  and  the  corresponding  values  of  are  given  for 
porH3  =  1.26  sec  (for  «r=  10-4  mhos/m,  H  =  90  km).  These  numerical  values 
for  a  =  2  are  considered  representative  of  the  daytime  near-equatorial 
ionosphere.  Time  histories  of  the  magnetic  field  in  the  cavity  are  shown  in 
Figure  3.  8  for  this  value  of  pcf^for  H/l^  =  0,  I  and  2.  It  is  seen  that  for 
H/ha  =  2,  a  substantial  decay  occurs  in  a  time  of  about  one  second. 


Hall- effect  plasma 


For  the  case  of  a  Hall  effect  plasma,  the  magnetic  field  varies  accord¬ 
ing  to  a  relationsnip  of  the  form  given  by  Equation  (3.  6),  with  the  w- functions 
being  replaced  by  W-functions  as  defined  by  Equation  (B.8.35). 


3.3.2  Frequency  Response 


The  cavity  frequency  response  to  a  steady  state  sinusoidal  input  of 
magnetic  flux,  for  the  case  of  a  plasma  slab  without  Hall  effects,  is  given 
by  Equation  (B.  8,  21)  as 

H(iw)  =  (1  -  [  k  H  tan  kh^.]-1)"1  (B.  8.  21) 
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where  the  wave  number  k  is  given  by  Equations  (B.  2.  28)  through  (B.  2.  30) 
for  complex,  perfect  and  ohmic  plasma  slabs.  For  a  plasma  slab  with  Hall 
effects.  Equation  (B.8.21)  above  is  replaced  by  a  four- component  tensor  with 
each  component  obeying  an  equation  similar  to  this  equation  {see  Eq.  (B.  8.  37)). 

3.  4  ESCAPE  FROM  A  VACUUM  CAVITY  INTO  AN 
EXPONENTIAL  PLASMA  REGION 

Section  B.  8.  4  of  Appendix  B  considers  the  problem  of  the  escape  of 

a  magnetic  pulse  from  a  vacuum  cavity  of  height  H,  which  is  bounded  below 

by  a  perfect  conductor  (where  E  =  0)  and  is  bounded  above  by  a  perfect  plasma 

whose  properties  simulate  the  nighttime  ionosphere  as  follows  (see  Fig.  3.  9). 

This  perfect  plasma  consists  of  a  lower  part  of  height  ha  ,  which  is  a  uniform 

perfect  plasma  with  phase  speed  VmQ.  Above  this  is  a  semi- infinite  horizontally 

stratified  perfect  plasma  region  whose  phase  speed  first  increases  exponentially 

with  increasing  altitude  (with  the  scale  height  2p)  and  then  eventually  approaches 

a  constant  uniform  condition  at  the  phase  speed  V  ,  where  V  >>-V  . 

mi  m  l  mo 

(See  Eq.  (B.  6. 1)). 

3.4.1  Step- Function  Response 

First  consider  the  step  function  response  of  the  cavity  to  a  unit  mag¬ 
netic  flux  applied  at  time  zero.  This  step  function  response  is  given  by 
Equation  (B.  8.  55)  as 

B  =  2  R[2  C  exp  (i«  t)]  (3. 14) 

v  n  m  r  n  ' 

where 

»  =H  V  / P  (3.15) 

n  n  mo  ' 

R[  ]  designates  the  real  part  of  the  function, 

and  the  C  's  and  |  's  are  complex  coefficients  of  H/ho  and  V  /V  (see 
Sec.  B.8.4).  Numerical  values  of  parameters  for  these  equations  were 
computed  for  VmQ/Vm  =  0.  05  and  a  =  H/(2p)  =  0.  968  (from  H  =  290  km  and 
p  =  150  km);  these  values  for  the  first  four  modes  are  shown  in  Table  II. 
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Fig.  3.9  Escape  of  a  Magnetic  Field  Through  An  Exponential 
Plasma  Region 
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TABLE  II 


VALUES  OF  CONSTANTS  FOR  EXPONENTIAL  PLASMA  PROBLEM 

(  V  /V  =  0.  05  ;  H/(2p)  =  290/300  =  0.  968] 

mo  mi 


*  Computed  for  =  0  ,  an  approximation  which  is  accurate 

enough  for  small  values  of  V  /V 

b  mo  m  i 


Assuming  2p/V 


m  r> 


=  300/310  =  0. 968  sec. 


4-» 


f 


*s 


A  time  history  of  the  magnetic  field  decay  is  shown  in  Figure  3. 10  for 

a=  %/(2p)  =  220  km/300  km  =  0.73.  Two  abscissa  scales  are  shown,  one 

for  a  generalized  time  scale  V  t/H  and  one  for  an  absolute  time  scale 

mo  - 

based  on  the  characteristic  time  (2p)/VmQ  =  300  km/310  km/sec  =  0.  97  sec. , 
which  is  representative  of  alow  latitude  nighttime  ionospheric  condition.  The 
magnetic  field  is  seen  to  initially  decay  exponentially  (Eq.  (3.1)  applies  here), 
ana  then  tends  to  oscillate  and  decay  more  slowly  at  the  lowest  characteristic 
frequency.  After  only  one  or  two  cycles  of  this  oscillation,  all  higher  modes 
die  out  and  the  subsequent  response  varies  simply  as 


B  =  0.  662  exp  (-  t/32)  cos  (0.  755  t) 


where 


T^Vmot/P 


(3.16) 


(3.17) 


The  characteristic  time  T  for  this  final  exponential  decay  is  seen  to  cor¬ 
respond  to  t  =  32  or 
f 


T  =  64p/V 


mo 


(3.18) 


It  is  interesting  to  note  that  while  the  magnetic  field  response  occurs 
predominately  m  the  lowest  frequency  mode,  the  same  is  not  true  for  the  rate 
of  change  of  the  magnetic  field  since  the  series  obtained  by  differentiating 
Equation  (3.14)  converges  much  more  slowly  than  the  series  in  Equation  (3.14). 
The  first  term  in  the  differentiated  series  remains  the  largest  term,  but  is  not 
much  larger  than  the  second  term  for  moderately  long  time  durations. 

3.4.2  Frequency  Response 

The  frequency  response  function  for  the  vacuum  cavity  for  a  steady- 
state  sinusoidal  application  of  magnetic  flux  is  given  approximately  by 
Equation  (B.  8.  50)  as 


I 

i 

I 


1 

L 


C 

c 

be 

[ 


>  ti 


J  u 


38 


’sponsc  To  a  Step- Function  Input 


H(iu>)  -  1  +  [  a|  tan  (a  +  b| )  -  l] 


(B.8.  50) 


_  1 


where 


1 

=  2u)p/V 

mo 

(B.  6.  3) 

a 

=  tan-x[  J1(^)/JQ(^)] 

(B.5.6) 

a 

=  &H/p 

(B.8.48) 

b 

=  i  (ha/p  +  iwe);  «  <  o 

(B.8.43) 

e 

=  V  /V  «  1 
mo  mi 

(B.8.44) 

This  variation  is  shown  in  Figure  3.11  as  a  function  of  the  generalized 
frequency  parameter  4  =  2<Jp/VrnQ  for  a  =  hg/ (2p)  =  220/300  =  0.  73  and  for 
a  =  0  .  (The  curve  for  a  =  0  corresponds  to  the  problem  of  Section  3.  2  for 
escape  into  a  semi- infinite  perfect  plasma. )  Also  shown  in  this  figure  is  an 
absolute  frequency  scale  for  the  characteristic  time  (2p)/Vmo  =  300  km  j 
310  km/sec,  representative  of  a  low  latitude  nighttime  ionospheric  condition. 

It  is  seen  from  Figure  3.11  that  there  is  a  conspicuous  near- resonant 
condition  at  the  generalized  frequency  4  =  “wP^mo  55  755  and  that  there 

are  a  sequence  of  other  higher  frequency  progressively  more  strongly  damped 
resonant  conditions.  The  first  four  characteristic  frequencies  are  listed  in 
Table  II  according  to  Equation  {B.  8.  47).  In  order  to  interpret  these  frequencies 
in  terms  of  nighttime  ionosphere  variations,  they  are  plotted  against  geomag¬ 
netic  latitude  in  Figure  3.12,  under  the  assumption  that  the  Alfven  speed 
varies  with  geomagnetic  latitude  in  the  same  manner  as  does  the  geomagnetic 
field  for  a  simple  dipole  representation,  i.  e.  ,  for 

V  =  2Y  AirPX  +  4  cos9  X  (3.19) 

a  ae 

where 

X  is  geomagnetic  latitude 

Vac  is  the  Alfven  speed  at  the  equator  (X  =  0) 
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Fig.  3.12  Characteristic  Frequencies  of  the  Nighttime  Ionosphere 


Both  generalized  and  absolute  frequency  scales  are  shown  in  this  figure  for 
=  220/300  and  110/300  ,  and  frequency  variations  are  shown  for  both 
Alfven  and  modified  Alfven  waves.  (No  curves  are  shown  for  the  Alfven  made 
at  low  latitudes  since  the  present  model  is  not  applicable  there.)  It  is  noted 
that  the  frequencies  increase  with  geomagnetic  latitude,  and  are  moderately 
larger  for  the  modified  Aifven  mode  at  intermediate  latitudes.  The  absolute 
frequencies  for  the  first  four  modes  are  seen  to  range  from  about  0.  1  to  2.  5  cps. 

3.  5  DESCENT  OF  A  HYDROMAGNETIC  WAVE  INTO  A  VACUUM  CAVITY 

Section  B.4  of  Appendix  B  considers  the  problem  where  a  semi¬ 
infinite  rectangular  hydromagnetic  wave  front,  descending  through  a  perfect 
hydromagnetic  plasma,  encounters  a  vacuum  cavity,  the  lower  edge  of  which 
is  bounded  by  a  perfect  conductor  iFig.  3.  13).  The  initial  descent  speed  of 
the  wave  front  depends  on  the  polarization  of  the  incide.nt  magnetic  field.  If 
its  magnetic  vector  is  in  the  plane  formed  by  the  ambient  magnetic  field 
vector  (Bq)  and  the  vertical,  the  wave  is  a  modified  Alfven  wave  traveling 
at  the  Alfven  speed  V  ;  if  its  magnetic  vector  is  perpendicular  to  this  plane, 

cl 

the  wave  is  an  Alfven  wave,  traveling  at  the  phase  speed  V  cos  <p ,  where  <p 

a 

is  the  acute  angle  between  the  ambient  magnetic  field  vector  and  the  vertical. 

In  either  case,  the  magnetic  field  produced  in  the  cavity  by  the 
descending  wave  increases  with  time  as  shown  in  Figure  3.  14,  according  to 
the  equation 

B  =  2B.[l  -  exp(-  V  t/H)]  (B.4.6) 

v  i  r  m 

where  Vm  is  given  by  Equation  (3.  2)  and  B.,  is  the  strength  of  the  incident 
magnetic  field. 

Descent  of  the  incident  wave  into  the  cavity  a.lso  produces  a  reflected 
wave,  whose  time  history  is  of  the  form  indicated  in  Figure  3.15.  It  may  be 
noted  that  the  reflected  wave  initially  has  a  negative  peak,  followed  by  an 
exponential  transition  to  an  eventual  reflection  coefficient  of  plus  unity. 

*  For  2p/V  =  300  km/300  km/sec. 

SC  / 

t  The  ratio  hs/(2p)  =  110/300  is  believed  to  be  more  realistic  for  the  nighttime 
ionosphere. 
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igure  3.  14  Growth  of  the  Magnetic  Field  in  the  Cavity 


3.  6  REFLECTION  FROM  EXPONENTIAL  BARRIERS 


Section  B.  5  and  B.  6  of  Appendix  B  consider  two  similar  problems, 
where  a  hydr ©magnetic  wave,  traveling  in  a  uniform  perfect  hydromagnetic 
plasma,  suddenly  encounters  a  region  of  exponentially  increasing  Alfven 
speed  (Figs.  3.16  and  3.17).  In  Section  B.  5,  the  Alfven  speed,  is  assumed  to 
increase  without  limit  (Fig.  3.16),  and  in  Section  B.6  is  assumed  to  have  a 
finite  upper  limit  (Fig.  3.17).  In  either  case,  the  incidence  of  a  semi¬ 
infinite  rectangular  wave  front  on  the  exponential  barrier  produces  a  reflected 
wave  of  the  form  shown  in  Figure  3,  18,  where 

\rmo  is  the  initial  phase  speed  (Fig.  3,  16) 

2p  is  a  scale  distance  associated  with  the 
phase  speed  variation  (Eq.  (B.  5.  3)) 

If  the  final  Alfven  speed  is  considered  unlimited  (Section  B.  5),  the  conspicuous 
peak  in  this  figure  is  infinite;  if  the  final  Alfven  speed  is  limited  (Section  B.  6), 
this  peak  becomes  finite. 

It  may  be  noted  in  Figure  3.18  that  the  polarity  of  the  reflected  wave  is 

reversed  from  that  of  the  incident  wave,  this  behavior  being  characteristic 

of  reflection  from  regions  of  greater  phase  speeds.  The  reflected  wave 

intensity  is  seen  to  increase  with  time  up  to  a  peak  value  at  the  time  4p/V  Q, 

after  which  time  the  wave  rapidly  approaches  the  endpoint  condition  of  a 

reflection  coefficient  of  minus  unity. 

Also  considered  in  Appendix  B  is  the  energy  transmitted  through  the 

exponential  barrier  for  the  case  of  incident  steady-state  sinusoidal  waves. 

The  fraction  of  the  incident  energy  transmitted  through  the  barrier  is  aero 

if  the  final  Alfven  speed  is  ’-nlimited;  is  equal  to  4V  /V  for  a  finite  phase 

mo  mi  r 

speed  limit  with  frequency  approaching  zero  (from  the  limiting  value  of  Eq. 

(B.  6.  5)),  and  increases  with  increasing  frequency  as  indicated  in  Figure  3.19. 


Fig.  3.  16  Unlimited  Exponential  Barrier  Problem 


Incident  wave 
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Fig.  3.  17  Limited  Exponential  Barrier  Problem 


3.  7  OTHER  PROBLEMS 


In  addition  to  the  material  specifically  described  above,  Section  B.  8.1 
of  Appendix  B  presents  some  additional  material  relevant  to  the  treatment 
of  more  complex  one-dirnensional  cavity  escape  problems  (Sec.  B.8.1),.  plus 
some  considerations  of  the  simple  problem  of  reflection  and  transmission 
of  hydromagnetic  waves  at  an  abrupt  barrier  betv/een  two  uniform  plasma 
regions. 


SECTION  IV 


APPLICATION  TO  THE  IONOSPHERE 


This  section  offers  some  brief  comments  on  the  applicability  of  the 
results  ~)f  the  preceding  sections  to  the  propagation  of  ultra- low- frequency 
(ULF)  waves  in  and  below  the  ionosphere-  Throughout  this  section  we 
disregard  the  fact  that  our  considerations  in  Section  II  were  limited  by  the 
assumption  of  an  isotropic  conductivity  law  for  the  plasma  regions  concerned 
and  are  thus  not  exactly  representative  of  the  anistropic  ionospheric  plasma. 
This  disregard  is  certainly  justified  for  the  essentially  qualitative  considera¬ 
tions  given  below;  the  reader  need  only  recognize  that  the  actual  ionosphere 
will  certainly  demonstrate  all  of  the  types  of  phenomena  characteristic  of  an 
isotropic  plasma,  plus  additional  complications  characteristic  of  the  aniso¬ 
tropic  conductivity  components. 

Consider  first  the  question  of  ULF  wave  propagation  via  the  sub¬ 
ionosphere  cavity,  such  as  may  be  excited  by  lightning  strokes  or  by  the 
magnetic  dipole  produced  by  the  detonation  of  a  nuclear  bomb.  Two  effects 
produced  by  a  magnetic  dipole  suddenly  appearing  in  or  below  the  ionosphere 
have  been  shown  in  Sections  LT  and  III.  First,  a  (hvdromagnetic)  wave  moves 
downward  from  the  dipole  origin  until  it  reaches  the  lower  edge  of  the 
ionosphere  and  then  immediately  (at  the  speed  of  light)  produces  a  world-wide 
magnetic  disturbance,  which,  in  turn,  then  relaxes  or  difuses  upward  into 
the  ionosphere.  This  relaxation  process  has  been  considered  in  detail  in 
Section  3.  4  for  a  fairly  realistic  model  of  the  nighttime  ionosphere  and  has 
been  shown  to  exhibit  various  resonant  frequencies  in  the  near  1  rps  region. 
Also,  the  relaxation  process  for  a  step  function  magnetic  flux  excitation  has 
been  shown  to  consist  of  one  dominant  frequency  below  one  cps,  which  persists 
for  tens  of  seconds,  plus  higher  frequency  overtones  which  decay  progressively 
more  rapidly.  r'or  the  daytime  ionosphere,  a  similar  but  more  complex 
response  would  be  expected,  added  complexity  being  introduced  by  the  Hall 
effect,  which  couples  the  Alfven  and  modified  Alfven  wave  components  in  the 
E-layer. 


It  may  be  inferred  from  the  above  disevssion,  that  the  distant  EJLF 
effects  of  lightning  strokes  will  also  be  subject  ;o  modification  in  the  same 
manner  as  for  nuciear  bursts  effects,  the  ouiy  basic  differences  for  these 
two  processes  being  in  the  frequency  spectrum  and  geographic  distribution 
and  polarization  of  the  primary  exciting  magnetic  field  traveling  through  the 
sub- ionosphere  cavity.  In  the  case  of  lightning,  of  course,  only  part  of  the 
wave  energy  occurs  within  the  ULF  range,  to  which  the  present  discussion 
is  limited. 

Regarding  wave  propagation  in  the  ionosphere  above  the  dense  (E-  layer) 
part  of  the  ionosphere,  it  has  been  shown  in  Sections  2.4  and  2.  5  that  the 
field  of  a  suddenly  created  magnetic  dipole,  when  incident  on  an  abrupt 
discontinuity  in  plasma  density,  can  create  a  surface  wave  of  strong  intensity 
traveling  along  the  interface.  Since  the  upper  ionosphere  contains  relatively 
strong  density  gradients  at  altitudes  of  about  400  kilometers  and  above,  it 
may  be  expected  that  such  a  form  of  wave  propagation  will  be  important  for 
the  upper  ionospheric  travel  of  bomb-produced  magnetic  disturbances.  It 
might  also  be  noted  that  this  surface- wave  type  of  wave  propagation  differs 
from  the  purely  radiation- type  propagation  associated  with  the  duct  theory 
of  upper  ionosphere  wave  propagation. 
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SECTION  V 


CONCLUDING  REMARKS 


This  report  has  presented  a  series  cf  solutions  of  one-  ar.d  three- 
dimensional  problems  concerning  electromagnetic  and  hydromagnetic  wave 
propagation  in  various  combinations  of  horizontally  stratified  plasma  and 
vacuum  regions,  chosen  to  simulate  significant  characteristics  of  the 
ionosphere  and  sub- ionosphere  cavity  at  ultra- low-frequencies  (near  1  cps). 

It  is  believed  that  the  results  of  these  studies  offer  an  increased  understanding 
of  propagati or  of  electromagnetic  and  hydromagnetic  wav  es  produced  by 
naclear  detonations  and  by  distant  lightning  strokes.  In  particular,  the 
frequency  response  and  step- function  response  of  the  nighttime  ionosphere 
to  ultra- lov. -frequency  waves  has  been  clarified  by  a  fairly  realistic  simula¬ 
tion  of  the  ionosphere.  The  corresponding  analysis  for  the  daytime  ionosphere 
would  be  considerably  more  difficult  (because  of  the  need  to  consider  Hall 
effects),  but  could  be  performed  in  a  similar  manner,  using  methods  outlined 
in  thir  report. 


APPENDIX  A 


THREE-DIMENSIONAL  ELECTROMAGNETIC  FIELD  PROBLEMS 


A.  1  INTRODUCTION 


This  part  of  this  report  is  concerned  with  the  calculation  of  the 
transient  electromagnetic  fields  produced  by  magnetic  dipoles  in  various 
types  of  horizontally  stratified  plasma  and  un-ionized  media.  It  is  assumed 
throughout  that  Maxwell's  equations  are  applicable  with  an  isotropic 
conductivity  law  and  that  the  displacement  current  may  be  neglected  in 
vacuum  regions.  This  assumption  of  isotropic  conductivity  is,  of  course, 
not  generally  applicable  to  hydromagnetic  wave  propagation.,  which  is  our 
principal  concern  in  this  report.  However,  to  a  considerable  extent,  both 
qualitatively  and  quantitatively,  the  electromagnetic  problems  considered 
herein  do  offer  some  insight  into  the  analogous  hydromagnetic  problems. 

The  following  electromagnetic  problems  are  considered,  after 
presentation  of  basic  equations  in  Section  A.  2.  The  magnetic  field  of 
a  magnetic  dipole  in  a  uniform  infinite  plasma  is  considered  in  Section  A.  3. 
Sections  A.  4  and  A.  5,  respectively,  consider  the  electromagnetic  fields 
produced  by  horizontal  and  vertical  magnetic  dipoles  located  in  a  semi¬ 
infinite  plasma  above  a  horizontal  interface  separating  the  plasma  above 
from  a  semi-infinite  vacuum  region  below.  Section  A.  6  considers  the 
electromagnetic  fields  produced  by  a  horizontal  magnetic  dipole  located 
in  a  semi-infinite  plasma  above  a  horizontal  vacuum  cavity,  the  lov/er  edge 
of  which  is  the  surface  of  a  perfect  conductor. 

For  each  of  the  above-mentioned  problems,  the  analysis  begins  with 
a  determination  of  the  electromagnetic  field  for  steady- state  sinusoidal 
variations  of  the  magnetic  dipole  moment.  This  is  followed  by  evaluation 
of  the  transient  magnetic  fields  produced  by  the  sudden  step-function 
creation  of  a  magnetic  dipole  moment  (see  Appendix  E).  Numerical  results 
are  presented  and  discussed  in  Section  II  of  the  text. 
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A.  2  BASIC  EQUATIONS 


A.  2.1  Maxwell’s  Equations 


The  basic  equations  used  in  this  Appendix  are  as  follows.  From 
Maxwell’s  equations  in  MKS  units 

VXB  =  |1 


where 


V  x  E  =  -dB/dt 
Jx  =  J  +  c~a  p.-1  oE/ot 


(A.  2.1) 
(A.  2. 2) 
(A.  2. 3) 


is  the  electric  field  vector 
is  the  magnetic  field  vector 
is  the  permeability  in  vacuum 
is  time 

is  the  velocity  of  light  in  vacuum 
is  the  current  density  vector 

is  the  total  current  density  vector,  including  the 
displacement  current,  which,  is,  however,  assumed 
zero  in  vacuum  regions 


We  consider  pJasma  regions  with  isotropic  conductivity  obeying  an  Ohm's 
law  relationship  for  sinusoidal  variations  of  the  form 


=  <r  E  ;  £  E  ~  exp  (ixot)j 


(A.  2. 4) 


where 


cr  is  the  conductivity 


Four  types  of  media  are  considered,  with  cr  defined  as  follows  for 


each. 


as  \c 

I  h 


m 


Perfect  plasma: 

tr  =  iu  /(fiV2) 

Complex  plasma: 

<r  =  co  f  itd/fpV2)  =  {  v  +  Lj)/(|iV2) 

Ohmic  plasma: 


where 


cr  =  cr 


cr  is  the  real  component  of  the  conductivity 

V  is  the  wave  propagation  speed  for  a  =  0 

v  is  a  characteristic  damping  frequency,  related 

to  cr  as 
o 


{A.  2. 5) 


(A.  2.6} 


(A.  2.7) 


v  =  a  V2i 


Vacuum 


cr  =  0 


(A.  2. 8) 


(A.  2. 9) 


The  electric  and  magnetic  fields  for  steady-state  sinusoidal  oscillations 
according  to  Equations  (A.  2.1)  through  (A.  2. 4)  may  be  expressed  in  terms  of 
a  vector  potential  F  as  (Wait  and  Campbell  1953  ) 


B  =  -  V2F  +  V  § 
E  =  -  iw  VXF 
§  =  V-F 
Y  =  ik 

k  =  \i-  icj  p.  cr 


(A.  2.10) 
(A.  2.11) 
(A.2.12) 
(A.  2.1  3) 

(A. 2.14) 


*  The  F  symbol  used  herein  differs  from  Wait  and  Campbell's  (1953)  symbol 
by  the  factor  iw  . 


59 


where 


is  a  scalar  potential 

is  the  wave  number 

is  the  propagation  constant 


and  where,  for  problems  with  symmetry  with  respect  to  the  x  =  0  plane  of 
an  XYZ-rectangular  coordinate  system  (Wait  and  Campbell,  1953) 


F  =  0 
x 


F  ,  Fz  ~exp  (La t  ±  u  z)  J^(X  p)  cos  n  0 


(A..  2.15) 


where 


=  vxTTy2 


(A.  2.16) 


J  is  a  Bessel  functicn  of  the  first  kind 
n 

n  is  a  constant 

X  is  a  parameter 


p  =  vxHy2 


(A.  2.17) 


0  =  tan  _1  (x  /y) 


(A.  2.18) 


For  each  problem  considered  here,  the  exciting  disturbance  is  taken 
as  a  horizontal  or  vertical  magnetic  dipole,  which  is  characterized  in  a 
uniform  plasma  by  the  vector  potential  (Wait  and  Campbell  1953) 


F  =  C  exp  (-  Yr)/r 


(A.  2.19) 


C  =  pmo/(4ir) 


(A.  2. 20) 


mQ  is  the  dipole  moment 

r  is  the  distance  between  the  dipole  and  any  point  (x,y,  z) 


1 


where  F  is  a  horizontal  vector  (F^)  *OT  a  horizontal  dipole  and  a  vertical 

vector  {F  )  for  a  vertical  dioole. 
z7 

A,  1.2  Boundary  Conditions 

Boundary  conditions  to  be  satisfied  at  the  boundary  between  plasma 
and  plasma -vacuum  regions  are  continuity  of  the  horizontal  electric  and 
magnetic  fields,  or  in  terms  of  the  vector  potential,  for  a  horizontal  plane 
boundary,  the  following  quantities  mui  be  continuous  across  the  boundary 
(Wait  and  Campbell,  1953): 


F  §,  v2F  ,  3F  / 3z  continuous 

z,  1  y  y 


(A.  2. 21) 


At  the  horizontal  surface  of  a  perfect  conductor,  the  boundary  condi¬ 
tion  to  be  satisfied  is  the  vanishing  of  the  horizontal  electric  field,  or  in  terms 
of  the  vector  potential 


F  =  3F  /dz  =  0  (A.  2.  22) 

z  y  v  ' 

A.  3  MAGNETIC  DIPOLE  IN  A  UNIFORM  PLASMA 

A. 3.1  Frequency  Response 

This  section  considers  the  electromagnetic  field  of  a  magnetic  dipole 
in  a  uniform  plasma.  The  dipole  is  assumed  to  be  located  at  the  origin  of  an 
XYZ-axis  system  (Fig.  A.l)  with  the  dipole  axis  parallel  to  the  Y-axis.  The 
electromagnetic  field  for  this  problem  lor  sinusoidal  oscillations  may  be 
expressed  in  terms  of  a  vector  potential  having  only  a  y-component  as 
(e.g..  Wait  and  Campbell,  1953) 

Fy  =  C  exp  (-  Yr)/r  (A.  3.1) 

where 

r  =  n/x  ijf  y2  +  z r  (A.  3. 2) 

y  =  ik  (A.  3.3) 
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The  magnetic  field  components  then  follow  from  Equations  (A,  2. 10) 
and  (A.  2.12)  as 


=  2  C  (1  +  Yr)  exp  (-  Yr)  cos  0  /r3 

(A. 3. 4) 

=  c£l  -f  Yr  -r  (Yr)2Jexp  (-Yr)  sin  0/r3 

(A.  3. 5) 

B.  =  0 

(A.  3.6) 

where  r,  0,  4*  ar-  spherical  polar  coordinates  with  respect  to  the  dipole  axis 
(Fig.  A.l).  These  three  equations,  taken  with  Equations  (A.  2.13}  and  (A.  2.5) 
through  (A.  2. 7)  of  Section  A.  2,  give  the  magnetic  disturbance  or  frequency 
response  in  a  plasma  resulting  from  a  steady-state  sinusoidal  variation  of 
magnetic  dipole  moment. 


A,  3.  2  Step-Function  Response 

It  is  now  desired  to  obtain  the  transient  response  of  the  plasma  to  a 
step  function  dipole  input.  This  may  be  done  conveniently  by  replacing 
to  by  -is  in  the  preceding  equations,  where  s  is  considered  tc  be  the 
transformed  variable  of  the  Laplace  transform  (Churchill,  1944),  by  then 
dividing  the  field  expressions  by  s,  and  then  inverting  the  resulting  Laplace 
transforms. 

Complex  plasma: 

For  the  general  case  of  a  complex  plasma,  the  step  function 
response  corresponding  to  Equations  (A.  3.4)  and  (A.  3. 5)  may  be  obtained 
with  the  aid  of  items  1  and  2  in  Table  E.l  of  Appendix  E  as 

Br(t)  =  2(C/r3)  [[,(T  -  1)  +  o2t  I}  (a  \l tz  -  1)  +  aIQ(a  \Tt^  -  1)J exp(-  qt) 


BQ(t)  =  (C/r3)^6 j(t  -D  +  fl  +  a  +  ja2)  6(t  -  1)]  exp  (-a) 

+  r^a2(T  +  a  -2  t/{t2  -  1))  lT(a  n/t*-!)  +  a  (1  +  ar/ (t2  -  l))lQ(a^Tl  -  l)^-exp(-  a-) 


exp  {-  a)  +  c2  (t)|S  {t  -  1)^  sin  8 


(A.  3.8) 


Aj(t)  =  f  IT  (a  \;t2  -  1)  exp  (-  ar)  «1t 


(A.  3. 9) 


(a  \r  t2  -  1)  =  Ij  (a  \!t*  -  1)  S(t  -  1)  /(a  n/t2  -  1) 


(A.  3.10) 


t  =  Vt/r 


(A.  3.11) 


a  =  2*  r  V  /V  =  |"Pl<r0rV 


(A.  3.12) 


Perfect  plasma; 

For  the  case  of  a  perfect  plasma  (v  =0),  the  magnetic  field  is  given  by 
the  simpler  equations 


Br(t)  =  2(C  /r3}  j\s(r  -  1)  +  6(r  •  1)J  cos  0 


Ohmic  plasma: 


(A.  3.13) 


B0(t)  =  (C/r3)  [s{t  -  1)  +  8(t  -  1)  +  Oj(t  -  1)J  sin  0  (A.  3.14) 


For  the  case  of  an  ohmic  plasma 


Br(t)  =  (2C/r3) j^erfc  Xj  +  xj  erf  Xjj 
BQ(t)  =  (C/r3)  £erfc  Xj  +  Xj  (1  +  2  x'j  )  erf  xj 


(A. 3.15) 


(A.  3. 16) 


i 

1 

I 


where 

x,  =  j  r  =  “(r/V)  4v7t 


(A. 3.17) 


I 


erf  Xj  =  (2/\J~rt)  exp  (-  xf) 


{A.  3.  18) 


I 


IS 


Transient  magnetic  field  variations  according  to  Equations  (A. 3. 7) 
through  (A.  3.18)  are  presented  and  discussed  in  Section  3.2  of  Section  IH: 

A.  4  HORIZONTAL  MAGNETIC  DIPOLE  IN  A  SEMI- INFINITE  PLASMA 

A. 4.1  Frequency  Response 

This  section  considers  the  electromagnetic  fields  produced  by  a 
horizontal  magnetic  dipole  located  above  a  horizontal  interface  separating 
a  plasma  (above)  from  a  vacuum  (Fig.  A.  2).  This  problem  has  been  analyzed 
by  Wait  and  Campbell  (1953)  for  sinusoidal  variations  of  the  dipole  strength, 
whose  results  for  the  magnetic  field  may  be  stated  in  the  form 


B  =  V  $  -  y2  F 


(A.  4.1) 


where,  for  the  plasma  region  (z>0) 


F  =  0 
x 


Fy  =  c[p(|z  -h|,p)  -  P(z  +  h,  p  )J 

Fz  =  -  2CV'2  £{83/9y  8z2)  N  (z  +  h,  p)  +  (92/9y  9z)  P  (z  +  h,  p)J 

$  =  co/ay)  jpflz  -  h  |  , p)  -  p (z  +  h, p) 

-  2  Y2  (92/9z2)  jjP(z  +  h,  p)  +  (a/az)  N  (z  +  h,  p/J  J 


HA.4.2) 


and  where  the  P  and  N  functions  are  defined  by  Equations  (C.  5)  and  (C.  7)  of 
Appendix  C,. 


1 

1 
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PlasmG  for  z  >0 
Vacuum  for  z<0 

Fig.  A.  2  Horizontal  Dipole  in  a  Semi-infinite  Plasma 

os  Designate  dipole  locations 
x  s  Designate  observation  points 


{ a.)  Reference  Problem  ( b.)  Reciprocal  Problem 


Fig.  A.  2  Two  Equivalent  Problems 
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It  may  be  noted  from  the  form  of  Equations  {A. 4.1)  and(A.4.2) 
that  the  following  reciprocity  law,  illustrated  in  Figure  A.  3,  is  valid  m  the 
plasma.  Specifically  if  z  has  the  same  numerical  value  as  h  and  either  h  or  z 
is  zero  the  magnetic  fields  for  these  two  cases  are  identical.  That  is,  the 
magnetic  field  at  the  plasma- vacuum  interface  created  by  a  dipole  at  the 
altitude  h  (Fig.  A.  3a)  is  the  same  as  the  magnetic  field  created  at  the 
altitude  h  by  a  dipole  at  the  interface  (Fig.  A.  3b). 

The  rest  of  this  section  is  restricted  to  the  special  case  where  the 
dipole  is  at  the  inferface  between  the  plasma  and  the  vacuum  (h  =  0).  For 
this  case  the  horizontal  magnetic  field  components  at  the  interface  are 
obtained  by  reduction  of  Equations  (A.  4. 2)  and  (A.  4.3)  as 


Bp  =  2  C( cos  0  /p3)j2  -  li.fi1  +  (f  +  5  +  12/f  4  12/fz  )  exp  (-  f)J 
Bg  =  C(sin  0/p3)  ^2  -  6/f2-f  (2  +  6/f  4  6/f  2)  exp  (-  f)J 


(A.  4.  3) 


where 


f  =  Yp  =  ikp 


(A.4.4) 


and  where  Bp  and  B^  are  the  radial  and  tangential  components  of  the  magnetic 
field  in  polar  coordinates  in  the  plane  where  z  =  0  (Fig,  A. 4.1).  These 
equaticrs,  taken  with  (A,  2. 5)  through  (A.  2.  8)  give  the  magnetic  field  frequency 
response  resulting  trom  a  steady-state  sinusoidal  variation  of  magnetic  dipole 
moment. 


A. 4.  2  Step-Function  Response 


It  is  next  desired  to  obtain  the  transient  magnetic  disturbance  in 
the  plasma  resulting  from  a  step  function  dipole  input.  This  transient 
response  function  for  perfect  and  ohmic  plasmas  may  be  obtained  from 
Equation  (A.  4.  3)  with  the  aid  of  Churchill’s  (1944)  Appendices  II  and  III  as 


B  =  2CB’'  cos  9  /p3 
P  P 

(A,  4.  5} 

Bg  =  CBg  sin  e/p3 
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where  and  are  normalized  values  of  the  magnetic  field  components 
(unity  at  t  -*•  co)  which  are  defined  below 


Perfect  plasma  (y  =  =  0) 


Bq  =  2  -  6t2  +  (6 t2  -  1)  S  (t  -  1)  +  6{t  -  1) 


Bg  =  2  -  3t2+  (3  t2  -  1)  S  (t  -  1) 


=  Vt/p 


(A.  4. 6) 


i 

1 

I 


H  9 


Ohmic  plasma: 


B*  —2  -  12  t j  erf  (|-  /VtJ)  -  erfc  (|-A/t7  ) 


+  (12\f  t  j/tt  +  1/n/tt  Tj  )  exp  (“~j  /tj) 


>  (A.  4.  7) 


i 


B0  =  2  -  &Tj  erf  (I/n/tj  )  -  erfc  +  6-vj  n/ir  exp  (-  j/tj) 


I 


where 


'i  =  t/(po-or2} 


(A.  4. 8) 


Magnetic  field  variations  according  to  Equations  (A,  4.  6)  and  (A.  4.  7)  are 
presented  and  discussed  in  Section  2.3  of  Section  II  . 
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A.  5  VERTICAL  MAGNETIC  DIPOLE  IN  A  SEMI-INFINITE  PLASMA 


A.  5.  1  Frequency  Response 


This  section  considers  the  electromagnetic  fields  produced  by  a  vertical 
magnetic  dipole  located  above  a  horizontal  interface  separating  a  plasma 
(above)  from  a  vacuum  (see  Fig.  A.  4).  This  is  the  same  problem  treated 
in  Section  A. 4,  except  for  the  orientation  of  the  dipole  axis.  This  problem 
has  been  treated  previously  by  various  authors  (e.g..  Wait,  1951;  Bhatta charyya 
1959,  1953),  who  have  established  basic  mathematical  approaches  to  the  prob¬ 
lem  and  have  dealt  in  detail  with  special  cases. 

The  electric  and  magnetic  fields  for  this  problem  may  be  expressed  in 
terms  of  a  single  magnetic  potential  F  as  (Bhatta charyya,  1963) 


B  =  GF 
z 


b  =a2F/aPaz 


E  =  ifcj9  F/9  p 
6 


(A.  5. 1) 


E  =  E  =  B  =  0 
p  z  9 


J 


where  G  represents  the  operator 


g  =  -  p-^a/apjpa/ap 


(A.  5.  2) 


and  where 


p  =  \fx2  +  y  2 


(A.  5.  3) 


9  is  the  azimuthal  angle  of  the  radius  vector 


and  where  the  subscripts  z,  p,  9  designate  the  corresponding  components  of 
the  electromagnetic  field  in  cylindrical  coordinates  (Fig.  A.  4). 

The  magnetic  vector  potential  F  has  been  obtained  in  an  integral  form 
by  Bhatta  charyya  (1963),  whose  Equations  (6)  and  (7)  for  the  upper  half-space 
(when  conduction  and  displacement  currents  are  neglected  in  the  lower  half¬ 
space)  may  be  reduced  to  the  explicit  form 


wnere 


Fzi  =  c|”p(jz|  ..p)  -5-  P(z1jP)  -5-  (2/y2)G  ^P(zj,  p) 

']] 


t  3N(z!,  p)  /  3zj 


for  z  +  h  >  0 


P(zi,p)  =  exp  (-Yrj)/r  j 

N{z„p)  =  IQ[  -£y  (r,  -  z,)J  Ko[  j;  Y  (r,+  z,)] 

rj  =  <s/p2  +  zf 


(A.  5. 4) 


(A.  5.  5) 


z1  =  z  +  2  h 


(A.  5. 6) 


Hereafter  we  restrict  our  attention  to  the  special  case  where  the 
dipole  is  at  the  interface  between  the  two  media  (h  =  0),  and  for  this  case, 
after  considerable  reduction,  Equations  (A.  5.1)  through  (A.  5.  6)  lead  to  the 
following  expressions  for  the  three  non-zero  components  of  the  electromagnetic 
field  at  the  boundary  between  the  two  media  (z  -  0). 


Bz(iu)  =  Bzs  Bz  (iw}  \ 

Bp  =  BZS  Bp  <«> 

E0(iw)  =  Bzs  Eq  (iw)  / 


(A.  5.  7) 


where 

Bzs  =  ”  C^p3  =  “  M-m0/(4Trp3) 


(A. 5. 8) 
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B*(L>)  =  18/(Y2p2)  -  [^/(Y^2)  +  18/ (Yp)  4  8  +  2  YpJ  exp(-  Yp)  (A.  5.  9) 


Bp(iw)  =  -  2p3  0/3p}p‘,(a/ap)[li4Yp)Ki4Yp)J 
E0  (i“)  =  2 iup|^3/{Yp)2  -  fl  +  3/{Yp)  +  3/{yp)2lexp  (-Yp) 


(A.  5. 10) 

(A.  5*  11) 


The  above  expressions  have  been  normalized  in  terms  of  the  vertical 
magnetic  field  for  a  static  dipole  (B  ),  the  normalized  variables  being 
designated  by  asterisks.  (Eqs.  (A.  5.  9)  and  (A.  5.11)  were  previously  obtained 
by  Wait  (1951)  and,  in  a  more  general  form,  by  Bhattacharyya  (1959).) 


A.  5.  2  Step  Function  Response 

It  is  next  desired  to  obtain  the  transient  electromagnetic  disturbances 
produced  by  a  step-function  dipole  input  for  a  perfect  plasma.  These  response 
functions  may  be  obtained  from  Equations  (A.  5.  9)  through  (A.  5.  11)  with  the  aid 
of  item  15  in  Table  E.  1  of  Appendix  E,  together  with  elementary  Laplace  trans¬ 
forms  (Churchill,  1944),  as 

B*  (r)  =  9  t2  -  (9t2  -  1)  S  (t  -  1)  -  2  6  (r  -  1) 
z 

B*(t)  =  2t£q  ^  (2  -r2  -  1}  -  (8  t2  -  7)Q  ^(2  T2  -  1)]  /[-nj-r2  -  l|]  (A.  5. 12) 

“  T  *2 

E*  (T)  =  2  V  £  3  r  -  6  (t  -  1)  -  3  T  S(t  -  1)] 


where 


T  =  Vt/p  (A.  5. 13) 

and  /here  the  Q’ s  are  Legendre  Emotions  of  the  second  kind,  which  may  be 
expressed  in  terms  of  tabulated  elliptic  integrals  by  the  relationships  given  by 
Abramowitz  et.al.'  (1964).  and  Erdelyi,  et.  al.  (1953,  I). 

Magnetic  field  variations  according  to  Equations  (A.  5. 12)  are 
presented  and  discussed  in  Section  2.4  of  Section  II  . 

❖  Equations  8.  3.  10  and  8.  3.  12  of  this  reference  appear  incorrect  and  were 
not  used. 
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A.  6  HORIZONTAL,  MAGNETIC  DIPOLE  ABOVE  A  VACUUM  CAVITY 


A.  6. 1  Frequency  Response 


This  section  considers  the  three-dimensional  problem  where  a 
horizontal  magnetic  dipole  is  placed  in  a  plasma  above  a  horizontal 
vacuum  cavity,  the  lower  edge  of  which  is  the  surface  of  a  perfect  con¬ 
ductor  (Fig.  A.  5). 

The  x,  y,  z-components  of  the  vector  potential  F  for  this  problem 
for  sinusoidal  oscillations  may  be  written  by  a  straightforward  generalisa¬ 
tion  of  Wait’s  (1953)  equations  as 


F 

yp 


F 

zp 


F 

zv 


F  =  0 
x 

00 

C  J  £exp(-  ujz  -  hj  )  +  fpexp(-uz)J  (X/u)  JQ(Xp)  dX 
o 

oo 

C/[’v+  exp  (+  Xz)  +  fv  exp  (-  Xz)J  JQ(Xp)  dX 
o 

00 

C  J  gp  exp  (-  uz)  X2  Jj  (Xp)  dX  cos  8 
o 

oo 

C  /  £gv+  exP  +  gv_  exP  ("  ^z)]  dx  cos  6 

o 


\ 

(A.  6.  1) 

> 


J 


where 


u  =\l\z  -  k2  =  X2  +  Y^ 
Y  ■=  ik 


(A.  6. 2) 


The  f's  and  g's  are  functions  of  X  to  be  determined,  the  dipole  axis  is 
parallel  to  the  Y-axis,  and 


v 

P 


is  a 
is  a 


subscript  designating  the 
subscript  designating  the 


vacuum  region 
plasma  region 
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location 

Origin  of  j 

Coordinates-^^ 


H 


Plasma  (z>0) 


Vacuum 


Fig.  A.  5  Horizontal  Dipole  Above  a  Vacuum  Cavity- 
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r  is  the  distance  from  the  dipole  to  a  point  P  (Fig.  A.  6.1) 
p  is  the  horizontal  projection  of  P 
0  is  the  angle  between  the  Y-axis  and  the  horizontal 
projection  of  the  line  OP  (Fig.  A. 5) 

Through  application  of  the  boundary  conditions  of  continuity  of  the 
electric  and  magnetic  fields  at  z  =  0  and  vanishing  of  the  horizontal  electric 
field  at  z  =  -H,  six  simultaneous  equations  are  obtained  for  the  unknown 
functions  fy±  ,  g v±  ,fp  and  g  .  Substitution  of  the  resulting  values  of  these 
functions  into  Equation  (A.  6.1)  yields  the  following  expressions  for  the 
vector  potential  and  its  divergence  §  at  z  =  0. 

00 

Fyp  =  C  /  [exP  (~  u|z  -  h  {  )  -  exp  (-u  (z  +  h))j  (X/u)  Jq  (Xp)  dX 
o 

CO 

Fyv  =  2  C  /ex P  ('  «b)  [cosh  (X(z  +  H))/sinh  Xh]  Jq  (Xp)  dX 
o 

oo  (A.  6.  3) 

Fzp  =  2  C/  [exP  ("  u(z  +  h))/(\  +  u  tanh  X  H)]  X  J,(Xp)  dX  cos  0 

o 

00 

zv  -  2  ^  J  |®XP  (-  uh)  smh  (X(z  +  H))/ [(x  +  u  tanh  XH)  sinh  XhJ  J  X  Jj  (Xp)  dX  cos  0 


=  2  C  V  cos  0 

z  =  0 

where 

co 

V  =  (d/dp)  J  v  Jq  (Xp)  dX 
o 

v  =[u/(X  +  u  tanh  XH)J  exp  (-  u  b) 


(A.  6. 4) 

(A.  6.  5) 

(A.  6. 6) 


$  =  V-  F 

P 
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It  is  now  desired  to  determine  the  horizontal  magnetic  fieldin  the 
vacuum  cavity.  This  will  be  done  first  by  considering  the  magnetic  field 
in  the  plasma  at  the  lower  edge  of  the  plasma  region  (z  =  0  ),  which  is 
given  in  polar  coordinates  in  the  plane  z  =  0  by  the  expressions 


Bp  =  a$/aP  (A.  6.  7) 

B0  =  p-1  d  /S9  (A.  6.8) 


where  as  given  by  (A.  6.4),  is  the  potential  of  the  horizontal  magnetic 

field.  By  combining  Equations  {A.  6,  4)  through  (A.  6.  8)  and  performing  the 

indicated  differentiations  with  respect  to  p  and  8  ,  the  following  expressions 

are  obtained  for  B  and  B~ 

P  0 


Bp  =  2Cp“3  cos  0  [Tx  -  T0]  (A.  6.  9) 

Bq  =  2Cp~3  sin8  [Tx  ]  (A.  6.10) 

where  °o 

T0  =  J  F  Xp  Jo<\p)  d  (Xp)  (A.  6. 11) 

o 

oo 

Tx  =  J  FJj  (\p)  d  '\p)  (A.  6.12) 

o 

F  =  Kpv  (A.  6. 13) 


In  order  to  evaluate  the  integrals  T0  and  Tx  ,  the  subsequent  analysis 
is  restricted  to  the  study  of  conditions  far  from  the  dipole  in  the  sense  that 
p>:>H  and  p  >>h.  For  these  conditions  the  function  F  behaves  as  follows 
(see  Eqs.  (A.  6.  6)  and  (A.  6.13)). 

r 


F  «  [vp/(l+  y  H)]  exp  (-  yh) 


\h<<  1 
\H<<  1 

\p<  L 


(A. 6. 14) 


f  * **  py  \2  +  v2 


;  X  «  k 


(A.  6.15) 


- 

«  p  Jl  k  (X  -  k) 

where  JL  is  some  large  number  (Lr  >>1).  For  values  of  X  not  covered  by  the 
above  expressions,  Xpv  is  a  slowly  varying  function  in  the  sense  that 
f'(x)/f(x)  «1. 

With  the  aid  of  the  above  observations,  the  integral  Tx  (Eq.  A.  6. 12: ) 
may  be  approximately  integrated  as  follows 

L 

Tx  «  J  [y  p/a  +  yH)]  exp  (-  yh)  Jx  (Xp)  dX 
o 

OO 

+  J  F  jf  (Xp)  d  (Xp)  (A.  6. 16) 

L 

j. 

v 

where  J,  (Xp)  is  an  asymptotic  expression  for  Jx  (Xp)  for  large  arguments 
(see  Eq,  (D.  18)  of  Appendix  D). 

The  first  term  on  the  right  hand  side  of  Equation  (A.  6. 16)  is  easily 
integrated,  the  second  term  may  be  integrated  approximately  by  Equation 
(1X20)  of  Appendix  D.  The  result  of  these  integrations  is 

Tj  »  -  [yp/(l  +  yH)]  exp  (-yh)  [J0(L)  -  l] 

+  [yp/(l  +  yH)]  exp  (-yh)  [Jo'CL)  +  (2L)-1  Jx  (L)] 

+  exp  (-yp)  (A.  6. 17) 

where  Equation  (A.  6.14)  was  used  to  evaluate  the  function  F  at  Xp  -  L.  The 
first  and  third  terms  in  Equation  (A.  6.17)  cancel, '  since  L  is  large,  so  that 

*  The  integration  procedure  used  here  to  evaluate  Tx  and  T0  requires  that 
kp  or  yp  >>1;  modifications  of  results  obtained  here  to  include  the  case  of 
smaller  kp  are  discussed  subsequently. 

_  i 

**  Their  sum  is  a  term  of  order  L  8  . 
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J0  (L,  ~  J0  (L);  the  fourth  term  is  negligible  compared  to  the  second  term. 
Hence  Equation  (A.  6.17)  reduces  to  the  form 


Ti  «[yp/(l  +  yH)]  exp  (-yh)  +  exp  (-yp)  (A.  6.18) 

Similarly  T0  is  found  to  be  given  by  the  expression 
T0»  Yp/{1+  YH)3  exp  {-  yh)  [L  Jj  (L)  -  0] 

-  [ Yp/(1  +  V  H)]  exp  (-y  h)  [  L  J*  (h)  +  £  J0"  (L)] 


-  yp  exp  (-yp)  (A.  6.19) 

The  first  and  third  terms  in  Equation  (A.  6.19)  cancel  '’  since  L  is  laige  and 
Ji  (L)»  Jx  (Lj.  The  fourth  term  is  of  order  L~^,  which  is  negligible  compared 
to  the  first  term  in  the  Tx  expression,  which  appears  beside  T0  in  the  only 
equation  in  which  T0  is  used  (Eq.  A.  6.  9).  Hence,  Equation  (A.  6.19)  reduces 
to  the  form 


To  »  -  yp  exp  (-yp)  (A.  6.  20) 

The  expressions  derived  above  for  T0  and  Tx  are  valid  for  all  large 
values  of  the  parameter  kp  but  not  for  smaller  values.  (The  used  expressions 
(D.  14)  and  (D.  15)  from  Appendix  D  require  that  the  parameter  A  be  a  slowly 
varying  function  near  x  =  xm  ,  v/hich  is  true  here  only  for  large  kp. )  Con¬ 
sequently,  another  set  of  expressions  must  be  derived  for  T0  and  Tx  for  smalle: 
values  of  kp  .  For  this  purpose,  it  is  convenient  to  approximate  Equation 
(A.  6.  6)  by  the  expression 

v  =  [u/X]  exp  (-uh)  (A.  6.  21) 


*  Their  sum  is  a  term  of  order  L 


n 

-t 
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which  is  valid  for  all  small  values  of  XH  (H<<:  p)  up  to  some  value  of  Xp  =  L 
where  L  is  a  large  number.  For  larger  values  of  Xp  ,  M)e  integrals  T0  and 
T:  (Eqs.  {A.  6. 11)  and  (A.  6.12))  do  not  depend  on  the  intermediate  shape  of  the 
integrand  function  F  {see  Appendix  D),  hence  Equation  (A.  6.  21)  can  be  used 
to  evaluate  Equations  {A.  6.11)  and  (A.  6.12)  for  the  entire  range  of  X  (for 
Yp  <  L  ;  yH  <<1).  The  result  is  {noting  Eqs.  (C.  5),  (C.  6),  and  {C.  9)  through 
(C.  11)  of  Appendix  C). 

co 

T0  =  p3  J Xu  exp  (-uh)  J0  (Xp)  dX 
o 


=  p3d2/db?  [(X/u)  exp  (-uh)  J0  (Xp)  dX 
=  p3  b2  P(h,  p)/bh2 


Ti  =  P2 


CO 

J  u  exp  (- 


uh)  (Xp)  dX 


f- 


-  p2d2/Sh?  /  u_1  exp  (-uh)  (Xp)  dX 


(A.  6.22) 


-  p2  o2  Q  (h,  p)/Shs 


(A.  6.23) 


Evaluation  of  the  derivatives  of  P  and  Q  and  neglect  of  higher  order  terms 
(h/p<<l)  gives  the  results 

T0  «  -  (1  +  yp  -  Y2  h2 )  exp  (-  Yp)  (A.  6.  24) 

Tx  ~  YP  exp  (-Yh)  +  (1  -  Y^2/p)  exp  (-Yp)  (A.  6.  25) 


We  now  have  two  sets  of  expressions  for  T0  and  Tj,  one  set  valid  for 
large  kp  (Eqs.  (A.  6. 18)  and  (A.  6.  20))  and  one  valid  for  kp  <,  L  ,  kH  <  <  1  and 


kh  <<  1  where  L  is  a  large  number  (Eqs.  (A.  6.  24)  and  A.  6.  25)).  Further¬ 
more  the  range  of  validity  of  the  expressions  overlaps,  hence  the  following 
solutions  valid  for  the  entire  range  of  X  may  be  obtained  by  inspection  of 
these  two  sets  of  equations 

To  ~  -  (1  -  y  p)  exp  (~yp)  (A.  6.  26) 

Tx  »  [yp/(l  +yH)l  exp  (-yh)  +  exp  (-yp)  (A.  6.  27) 


which  are  seen  to  satisfy  Equations  (A.  6.18)  and  (A.  6.  20)  for  large  yp  and 
Equations  (A.  6.  24)  and  (A.  6.  25)  for  kh  and  kH  small  (kp  not  necessarily 
small). 

The  horizontal  magnetic  field  components  at  the  upper  edge  of  the 
vacuum  cavity  then  follow  from  Equations  (A.  6.  9),  (A.  6.10)  (A.  6.  26),  and 
(A.  6.  27)  as 


B0  =  2  C  (sin  8  /P3)  Typ  exp  {-  Yh)/(1  +  YH)  +  g  exp  (-Yp)] 
Bp  =  4  C  (cos  0/p3)[-|  Yp  exp  (-  Yh)/(1  +  YH) 

+  (1  +  -jYp)  g  exp  (  -  Yp)J 


V  (A.  6.  28a) 


where 

g  =  1  (for  z  =  0) 


(A. 6. 28b) 


In  a  similar  manner,  the  magnetic  field  components  may  be  evaluated 
at  the  surface  of  the  perfect  conductor  (z  =  -  H),  the  only  essential  change  in 
the  derivation  being  that  the  potential  of  the  horizontal  magnetic  field  is  given 
by  expressions  like  Equations  (A.  6.  4)  and  A.  6.  5)  with  v  replaced  by 
v/ cosh  k  h  =  vl cos  yh.  The  resulting  expressions  for  the  horizontal  magnetic 
field  compon.  its  are  found  to  satisfy  Equations  (A.  6.  28)  provided  that  g  is 
defined  as 


g  -  sec  y  H  (for  z  =  -  H) 


(A. 6. 28c) 


80 


.f iVi-;p~  r  ... 


A.  6.  2  Step- Function  Response 

The  transient  response  of  the  vacuum  cavity  to  a  step  function  dipole 
input  may  be  obtained  by  applying  items  1,  2,  8,  23  and  .4  of  Table  E.  1  to 
Equation  (A.  6.  28)  and  i~  given  by  the  following  equations  for  the  three  types 
of  plasma  considered. 

BQ  =  2  C  (sin  0/p?)  |jp/H)  Bj  +  B2J  ^ 

Bp  =  4  C  (cos  8/p3)£^-{p/H)3j  +  B2  +  B3J  J 


Perfect  plasma: 


B,  =  exp  T-  (Vt  -  h)/Hj  S(Vt  - 


h) 


For  z  =  0: 


B2  =  S  (Vt  -  p) 


B3  =  6  (Vt/p  -  1) 


For  z  =  -H: 


B3  =  i  (p/H)  / cosh  -I  (Vt  -  p)/H 


Complex  plasma: 


Bj  =  exp  (-  ihv/V)  S  (Vt  -  h)  ;  Vt  -  h  <  0+ 


Ohmic  plasma: 


Bj  =  exp  (h/H  +  r  j)  erfc  £  'n/t1  +  \  (h/H)A/Tj  J 


B2  =  erfc  ( n/t"1  )  (for  z  =  0) 

B3  =  |  erf^  Wt-1  )/*Jtz  (for  z  =  0) 


(A.  6.  29} 


(A.  6.  30} 


(A.  6.  31) 


(A.  6.  32) 


B2  =  y  +  -rr  -1  sgn  (Vt  -  p)  tan-1  £sinh  |-|  it  (Vt  -  p) /H j  J  (A.  6.  33) 


(A.  6.  34) 


(A.  6.  35) 


(A.  6.  36) 
(A.  6.  37) 
(A.  6.  38) 
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(A.  6. 29) 


Ti  -  t/{p-o-0H2) 

~z  =  t/(^u  trop2)  (A.  6 . 40) 

Magnetic  field  variations  according  to  the  above  equations  are  presented 
and  discussed  in  Section  2.  5  of  Section  II  . 


APPENDIX  B 


- -  | 

! 

ONE -DIMENSIONAL  HYDROMA GNETIC  WAVE  PROBLEMS  j 

B.  1  INTRODUCTION 

This  appendix  is  concerned  with  the  analysis  of  four  one- dimensional 
problems  of  hydromagnetic  wave  propagation,  which  are  related  to  the 
electromagnetic  phenomena  produced  by  a  high  altitude  nuclear  detonation. 

The  following  problems  are  considered,  using  the  basic  assumptions  and 
equations  given  in  Section  B.2. 

Section  B.3  considers  the  problem  where  a  magnetic  disturbance 
is  suddenly  applied  inside  a  vacuum  cavity,  located  between  a  perfect 
conductor  and  a  semi-infinite  hydromagnetic  plasma.  The  decay  of  this 
disturbance  in  the  cavity  and  its  propagation  into  the  plasma  are  examined. 

Section  B.4  considers  the  problem  where  a  hydromagnetic  wave, 
descending  through  a  semi-infinite  plasma,  encounters  a  horizontal  vacuum 
cavity,  the  lower  edge  of  which  is  bounded  by  a  perfect  conductoi .  The  growth 
of  the  magnetic  field  in  the  cavity  and  the  strength  of  the  wave  reflected 
back  into  the  plasma  are  examined. 

Sections  B.5  and  B.6  consider  two  similar  problems  where  a 
hydromagnetic  wave,  traveling  in  a  uniform  plasma,  suddenly  encounters 
a  region  of  exponentially  increasing  Alfven  speed.  For  the  problem  in 
Section  B.5,  the  Alfven  speed  is  assumed  to  increase  indefinitely  with  increas¬ 
ing  altitude,  but  in  Section  B.  6  is  assumed  to  have  a  finite  upper  limit. 

Reflection  and  transmission  coefficients  are  calculated  for  these  problems. 

Section  B.7  considers  the  reflection  and  transmission  coefficients 
for  a  steady-state  sinusoidal  wave  crossing  an  abrupt  boundary  between  two 
different  plasma  regions. 

Section  B.8  considers  the  problem  where  a  magnetic  disturbance  is 
suddenly  or  sinusoidally  applied  inside  a  horizontal  vacuum  cavity,  located 
above  a  perfect  conductor  and  below  various  combinations  of  horizontally 
stratified  plasma  regions  simulating  the  lower  ionosphere.  The  decay  of  the 
magnetic  field  in  the  cavity  and  the  cavity  frequency  response  are  examined. 


B.2  BASIC  EQUATIONS 


The  basic  equations  used  are  as  follows.  From  Maxwell's  equations, 
with  the  displacement  current  being  neglected: 

Vx  E  =  -  9B  /3t  (B.2.1) 

VXB=iiJ  (B.2. 2) 

The  following  different  types  of  media  are  considered,  with  the  indicated 
relationships  between  the  electric  and  current  fields. 

Vacuum  region 

1=0  (B.2. 3) 

Simple  conductor  (hereafter  called  ohmic  plasma) 


J  =  trE 

Perfect  hydromagnetic  plasma 
E  =  -  v  XB 

o 

p 3v  /at  =  JXB 
r  o 


(B.2. 4) 


(B.2.  5) 
(B.2. 6) 


where 

v  is  the  plasma  velocity 

Bq  is  a  constant  ambient  magnetic  field  vector 


Complex  hydromagnetic  plasma 
E  =  -  v  XB 

o 


where 


p3v/  at  +  pvv  =  JXBq 
v  is  a  collision  frequency 


(B.2. 7) 
(B.2. 8) 


The  above  current-electric-field  relationships  are  standard  relation¬ 
ships  except  possibly  for  the  last  set.  These  particular  relationships,  for  a 
complex  hydromagnetic  plasma  are  applicable,  for  example,  for  hydromagnetic 
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wave  propagation  in  a  weakly  ionized  plasma  without  substantial  Hall  effects; 
in  particular  for  the  F-layer  of  the  ionosphere,  if  v  is  taken  equal  to  the 
frequency  of  collisions  of  positive  ions  with  neutral  particles. 


J  = 

0*  E  -  a- s  E  xbQ  -  {<tq  -  cr1}(E*bo)  bQ 

(B.2.9) 

cr 

o 

=  oo 

(B.2.10) 

b 

=  B  /  IbJ 

(B.2.11) 

o 

or 

o 

is  the  parallel  conductivity 

«i 

is  the  Pedersen  conductivity 

*3 

is  the  Hall  conductivity 

and  where  and  <r2  are  considered  to  be  real  constants.  These  equations 
and  assumptions  are  appropriate,  for  example,  to  the  lower  E -layer  of  the 
ionosphere,  above  about  90  km,  for  ^quencies  considerably  lov/er  than  the 
ion  collision  frequency. 

One -dimensional  condition 

This  appendix  is  restructed  to  the  treatment  of  one -dimensional 
problems  (3/3x  =  3/3y  =  0),  for  which  case  the  preceding  equations  take 
the  following  simplified  forms,  when  expressed  in  an  X,  Y,  Z  rectangular 
coordinate  system.  Equations  referring  to  the  modified  Alfven  mode  of 
wave  propagation  are  placed  below  on  the  left  hand  side  of  the  page,  those 
for  the  Alfven  mode  on  the  right  hand  side,  and  common  or  miscellaneous 
relationships  in  the  center.  It  is  also  assumed  here,  without  loss  of 
generality,  that  the  constant  ambient  magnetic  field  vector  (Bq)  lies  in  the 
YZ  -plane  (Bqx  =  0). 


Maxwell’s  equations 


Modified  Alfven  Mode 


Alfven  Mode 


3B  /at  =  -  as  / 3z 
y  x 


/at  =  o 

z 


8B  /at  =  +  3E  / 8z  (B.2.12) 
x  y 


uJ  =  -  3B  /  3z 
r  x  y 


J  -  0 
z 


=  +  3Bx/3z  (B.2.13) 


Complex  hydromagnetic  plasma 


E  =  -  v  B  +vB 
x  y  oz  z  oy 


E  =  ■*■  v  B  ;  E  =  -v  B  (B.2.14) 
y  x  oz’  z  x  oy  '  ' 


p(3v  /at  +  vv  )  =  -  J  B 

r'  y  y'  X  OZ 


p{3vz/at  +  wz)  =  +  JxBoy 


p(3vx/at  +  wx)  =  +  JyBoz  (B.2.15) 


(For  a  perfect  hydromagnetic  plasma  set  v  =  0  in  the  above  equations.) 


Ohmic  plasma 


J  =  a  E 
x  x 


J  =  c  E 

y  y 


(B.2.16) 


Hall -effect  plasma 


rx  =  n  ex  +  Er 


j.  =  -  E  +  O'!  E 


J  =  0 

z 

E  =  -  E  cos  0 

y  r 


(B.2.17) 


E  =  E  sin  0 
z  r 


where 


Jr  =  -  Jy  cos  0  (B.2.18) 

9  is  the  angular  deviation  of  the  ambient  magnetic  field 
vector  (Bq)  from  the  vertical  Z-axis  toward  the  Y-axis. 

is  the  resultant  electric  field  in  the  YZ -plane. 

Differential  equations 

The  differential  equations  for  B  and  E  obtained  from  the  above 
relationships  for  a  complex  or  perfect  plasma  are 


92E/9t?  +  v  9E/9t  =  V2  92  E/9z2 


m 


92B/9t2  +  v  92B/9t  =  V2  92  B/9z2 

m 


^complex  plasma)  (B.2. 19} 


where 


Modified  Alfven  mode 


Alfven  mode 


V  =  Jb* 
m  o 


/Urp) 


V  =  Jbz  /(pp) 
m  oz  'rr' 


(B.  2.20) 


and  V-  0  for  a  perfect  plasma.  Also,  for  vacuum  regions 
9B/9z  =  0  (vacuum) 
and  for  a  simple  ohmic  plasma 


(B.  2.21) 


92B/9zz  =  per  9B/9t 


92E/9z2  =p.or9E/9i 


ohmic  plasma  (B.2.22) 
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For  a  Hall-effect  plasma,  the  differential  equations  resulting  from 
the  coupled  Equations  (B.2.12),  (B.2.13),  (B.2.17),  and  (B.2.18)  maybe 
conveniently  expressed  in  terms  of  two  separated  mode  equations  (designated 
by  =  superscripts)  as  follows. 


where 


a2Ei  /3z2  =  uo-  *  dE±  /at 


as*  /at  =  aE~  /az 


B±  =  B  +  b*  B 
x  y 

Ei  =  E  -  b±E 
y  x 


(B.  2.23) 


(B.2.24) 


c r  =  c j  /cos2  0  +  <r2  b^  /  cos  0 


b~  =  3  /  (-  1  ±^1  -  &^) 

3  =  2  (<r2  /o-j^)  cos  0/sin2  0 


(B.2.25) 


and  where  the  inverse  relationships  to  those  of  Equation  (B.2.24)  are 


Bx  =  (b+  B"  "  b~  B+)/<b+  "  b~) 
By  =  (B+  -B“)/(b+-b") 

Bv  =  -  (£+  -  Z")/(b+  -b") 


E  =  (b+  E"  -  b"  E+)/(b+  -  b") 


(B.2.26) 


Sinusoidal  wave  motion 


For  sinusoidal  wave  motions  in  uniform  plasma  regions,  solutions  of 
Equations  (B-?.I9),  i-B.2.22),  and  (B.2.23)  take  the  form 


where 


B,  E  ~exp  i(d)tikz)  =  exp  (i  u>t  ±  y  z) 


k  is  the  wave  number 
y  is  the  wave  propagation  constant  (ik) 


(B.2.27) 
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Substitution  of  Equation  (B.2.27/  into  Equations  (B.2.19),  (B.2.22),  and 
(B.2.23)  gives  k  and  y  as 


r 


+  (i  «J)  V  /  V 

m 

complex  plasma 

(B.2,28) 

i«J  /  V 

m 

perfect  plasma 

(B.2.29) 

Vi  IJU 1-1  <T 

ohmic  plasma 

(B.2.30) 

y i  toM  <r~ 

Hall -effect  plasma 

(B.2.31) 

B,  3  ESCAPE  OF  A  MAGNETIC  PULSE  FROM  A  VACUUM  CAVITY 

B.3.1  Derivation  of  Equations 


This  section  considers  the  escape  of  a  finite  magnetic  pulse,  suddenly 
applied  to  <x  vacuum  cavity,  which  is  bounded  below  by  a  perfect  conductor 
and  above  by  a  plasma  medium  (Fig.  B.l). 

The  initial  conditions  of  this  problem  are 


where 


B  =  E  =  0  ;  t<  0  ;  z  >  0 


B(z)  =  Bvq  =  §q/H  ;  -  H  <  z  <  0  ;  t  =  0 


$o  =  ^  B(z,  0)  dz  =  BvqH 
-H 


(B.3.1) 


(B.3.2) 


§  is  the  initial  magnetic  flux  in  the  vacuum  cavity 
Bvq  is  the  initial  (uniform)  magnetic  field  in  the  cavity 
H  is  the  height  of  the  cavity 

z  is  distance  above  the  upper  edge  of  the  cavity 
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The  boundary  conditions  to  be  satisfied  for  t  >  0  are  continuity  of 
B  and  E  at  the  boundary  at  z  =  0,  vanishing  of  the  electric  field  at  the  sur¬ 
face  of  the  perfect  conductor  (z  =  -H),  and  vanishing  of  the  field  at  infinity 
(z  =  +  oo). 

This  problem  may  be  conveniently  solved  by  taking  the  Laplace 
transforms  of  the  appropriate  differential  equations  with  respect  to  time 
and  satisfying  the  boundary  conditions  for  the  transformed  variables. 

B.3,2  Complex  Hydromagnetic  Plasma 


Considering  first  the  case  of  a  complex  hydromagnetic  plasma, 
the  Laplace  transform  of  the  differential  Equation  (B.2. 19'  is 

s2E"(z,s)  +  v  sE'"(z,  s)  =  Vjna2E'r(z,s)/3z2  ;  z  >0  (B.3.3) 

where  the  component  subscripts  given  in  Section  B.2  are  omitted  here, 
asterisks  are  used  to  designate  Laplace  transforms  with  respect  to  t,  and  s 
is  the  transformed  variable  of  the  Laplace  transform.  The  solution  of  this 
equation  which  vanishes  at  positive  infinity  is  of  the  form 

E  (z,  s)  =  A(s)  exp  £  -  's/s2  +  v  s  (z. / ;  z  >0  (B.  3.4) 

where  A{s)  is  to  be  determined.  The  corresponding  expression  for  the 
magnetic  field  is  (using  Eq.  (B.2.  12)) 


B  (z,s)=±A(s)  si  sz  +  v  s  /(sVm)J  exp  j  -s!  sHvs  (z/V^lJ  ; 


z  >  0 


(B.3.5) 


where  upper  signs  apply  to  the  modified  Alfven  mode  and  lower  signs  to  the 
Alfven  mode. 

From  Equations  (B.  2.12)  and  (B.  2.21 ),  with  the  electric  field  vanishing 
at  z  =  -H,  the  Laplace  transform  of  the  electric  field  expression  for  the 
vacuum  cavity  takes  the  form 
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E%,  s)  =T(z  f  H)  [s  B*(s)  -  Byo  ]  ;  -  H<  z  <  0 


(B.3.6) 


where  B  (t)  is  the  magnetic  field  in  the  cavity  (independent  of  z). 

The  unknown  functions  A(s)  and  Bv(s)  are  determined  by  the 
requirement  of  continuity  of  E  and  B  at  z  =  0  as 


A(s)  r.  ±  V  HB  /(V  +  H  n/s4  x  vs)  (B.3.7) 

m  vo  m 

B*(s)  =  (Bvq/s)  [fL/s2  +  vs  /(Vm  +  H  «v/  s2  +  vs)j  (B.  3. 8) 


These  Laplace  transforms  may  be  placed  in  a  convenient  form  for 
inversion  by  rationalizing  their  denominators,  dividing  remaining  irrational 
terms  by  \js2  +  v  s  ,  and  applying  partial  fraction  expansions  to  rational 
polynomials  in  the  denominators.  For  example.  Equation  (B.3.8)  maybe 
rewritten  in  the  form 


B^(s)/Bvo  =  s'"1  +  (s-1  -  a  \!  (as  +  2  bj  /as  )  £(as  +  b  -  \lbz TT) 


(as  +  b  +  \l bz  +  1)  1  J  /(  2 -s/FTl) 


(B.3.9) 


where 


a  =  H/V 


m 


b  =  f  v  H/V 
*  m 


(B.  3.10) 


Inversion  of  the  above  Laplace  transforms  may  be  accomplished  with 
the  aid  of  standard  elementary  Laplace  transform  tables  (Churchill,  1944) 
and  Erdelyi’s  (1954,  I)  transform  5.3  (27).  The  resulting  expression  for  the 
magnetic  field  in  the  cavity  is 


B(t)/Bvo  =  exp(-  bt/a)  exp((3t/a)/(p  -  b)  +  exp(  •  pt/a)/((3  +  b) 


+  2  sinh(pt/a)  IQ(-  bt/a)  -  I+  -  I_J  /(2p 


(B.3  11) 
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where 


t/a 

I.  =  (0±b)  V  exp{“  pT)  I  (-  br)  dT  (B.3.12) 

o 

p  -  -JITF"  (B.  3.13) 

and  a  and  b  are  given  by  Equation  (B.3.10). 

B.3.3  Perfect  Hydromagnetic  Plasma 

For  a  perfect  hydromagnetic  plasma,  the  Laplace  transform  of 
the  magnetic  field  is  simply 

B%,s)  =|j3voH/(Vm  +Hsj]  exp  (-  sz/Vml  ;  z  >0  (B.3.14) 

the  inverse  transform  of  which  is  (Churchill,  1944) 

B(z,t)  =  Bvoexp£-  (t  -  z/H)J  S(t  -  z/H);  (B.3.15) 

where 

T  =  Vmt/H  (B.  3.16) 

and  where  the  magnetic  field  in  the  cavity  is  given  by  (B.3.15)  with  z  =  0. 

B.3.4  Ohmic  Plasma 

For  a  simple  ohmic  plasma,  the  Laplace  transform  of  the  magnetic 
field  expression  is  easily  found  to  be 

B%,  s)  =  ByoT  exp[-VsT"(z/H)]y/[VsT(l  +  \TsT)]  ;  z>  0  (B.3.17) 

where 

T  =  po-  H2  (B.  3.18) 
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and  the  inverse  o£  this  is  (Churchill,  1944) 


B(z,t)  =  Bvoexp{Tj  +  z/H)  erfc 


[n/t^  +  (z/H)/*/ti];  z  >  0 


where 


Tj  =  t/T 


(B. 3.19) 


(B.  3.20) 


The  magnetic  field  in  the  cavity  is  obtained  from  n.quation  (B.3. 19) 
with  z  =  0  as 

B  (t)  =  B  w(i7t7T)  (B.  3.21) 

V  YO 

where 

w{x)  =  exp  (-x2)  eric  (-  ix)  (B.  3.22) 

and  where  w(x)  is  a  tabulated  function  (e.g.,  Abramowitz,  et  al,  1964). 

B.3.5  Hall -Effect  Plasma 


For  the  case  of  a  Hall-effect  plasma,  the  Laplace  transforms  of  the 
differential  equations  (B.2.23)  with  respect  to  time,  for  the  plasma  region, 
in  terms  of  the  variables  E~  (see  Eq,  (B.2.24)),  are 

d2E=  (z,  s)/d  z2  =  jicr  ~  s  E~  (z,  s)  ;  z  >0  (B.3.23) 

The  solutions  of  these  equations  which  vanish  at  positive  infinity  are  of  the 
form 

i  i  / 

E  (z,s)=A  (s)  exp  (-N/lio--  s  z)  ;  z  >0  (B.3.24) 

where  the  A' s  are  to  be  determined.  The  corresponding  expressions  for  the 
magnetic  field  are  (using  Eq.  (B.2.12)) 

B"  (z,  s)  =  -  £  A*  (s)  /s  jexp  (-  J[l  a*  s  z)  :  z  >0  (B.  3. 25) 


- 


94 


The  corresponding  expressions  for  the  electric  field  *n  the  vacuum 
cavity  follow  from  Equation  (B.3.6)  as 


e'  (z,s)  =  -(2tH)[sBjv(s!-Byvoj  ^ 


>  -  H<z  <0 


(B.  3.26) 


where 


B 

xvo 


B 


yvo 


are  the  initial  XY  magnetic  field  components  in  the 
vacuum  cavity 


or  in  terms  of  E 


(see  Eq.  (B.2.24)) 


E±%,s)  =  (z  +  PI)  [s  ET*{s)  -  B^j 


where 


B“  =  B  +  b~  B 
vo  xo  yo 


(B.  3.27) 

(B.3.28) 


•  * 

The  unknown  functions  A=  (s)  and  B^  (s)  are  determined  by  the  re¬ 
quirement  of  continuity  of  E~  and  B~  at  z  =  0  (using  Eqs.  (B,  3.  24),  (B.  3. 25), 
and  (B.  3, 27))  as 


B“  (s)  =  B~  "  (0,  s)  =  H  B^q  j  H  s  +V s/m c~ 


(B.3.29) 


=  -  B~  "'  (s)  «/s/|iir* 


A"  (s)  =  -  B;  (s) 


(B.3,30) 


The  Laplace  transform  of  the  magnetic  field  in  the  plasma  then  follows  from 
Equations  (B.3.25),  (B.3.29),  and  (B.3.30)  as 

B~  (z,  s)  =  B~vo  T ,  exp  \l  s  T^  (z/H)J  /  \i  s  T^  ( 1  +  \l  s  T J  (B.3.31) 

where 

Ta  =  Mo-±H2  (B.3.32) 


%  Jg'.  j.  ~***-‘U**jt*P'  1 
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The  magnetic  field  in  the  plasma  is  obtained  by  inverting  Equation  (B,  3.31), 
the  result  being  the  same  as  Equation  (B.3.19),  with  ±  subscripts  or  super¬ 
scripts  applied  to  B  and  T.  For  the  magnetic  field  in  the  cavity  (z  =  0),  the 
result  is  (see  Eq.  (B.3.21)) 


B=  (t)  =  B“  w  (i  *n/  t/T ,  ) 
v  vo  = 


(B.3.33) 


or  in  terms  of  B  and  B  (see  Eqs.  (B.3.25),  (B.2.26),  and  (B.3.28)). 
x  y 


where 


Sx(t)  =  Cu  w+  +  w_ 


B  (t)  =  ca  'v+  +  C22  w. 


C,  =  -  c“(B  +  b^B  ) 

'  vn  vrr 


a2  =  c+(B  +  b"B  ) 
^2  '  xo  yo' 


(B.  3.34) 


Ca  =  (Bxo  +  i>+By0)/(b+  -  b-) 


C22  =  -  (B  +  b"B  )/(b  -  b") 
'  xo  yo 


c4  =  b4  /(b+  -  b") 


w  =  t/T^  ) 


•...  V  :  S!i  ;  -.A- X A-T- 


.  v 


B.  4  DESCENT  OF  A  HYDROMAGNETIC  WAVE  INTO  A  VACUUM  CAVITY 


This  section  considers  the  descent  of  a  hydromagnetic  wave  front 
from  a  perfect  hydromagnetic  plasma  into  a  vacuum  cavity  bounded  at  its 
lower  edge  by  a  perfectly  conducting  plane  (Fig,,  B.2). 

The  initial  conditions  for  this  problem  are 

B  =  E  =  0  ;  t<  0  ;  z  <  0  (B.4.1) 


B(z)  =  B_.(z) 

X. 

E(z)  =  0 


t  =  0  ;  z  >0 


(B.4.2) 


The  boundary  conditions  and  differential  equations  are  the  same  as  in  the 

preceding  section.  The  solution  to  this  problem  may  be  considered  to  consist 

of  three  parts:  an  incident  wave  of  amplitude  B^,  the  disturbance  in  the 

vacuum  cavity  of  amplitude  Bv>  and  a  reflected  wave  of  amplitude  B^,  whose 

respective  amplitudes  may  be  conveniently  evaluated  by  considering  the 

conservation  of  magnetic  flux  flowing  into  the  cavity.  The  net  inflowing 

flux  is  V  (B.  -  B  )  which  must  be  equal  to  the  rate  of  increase  of  flux  in 
mi  r  ^ 

the  cavity  HSB^/St  ,  or 

H9Bv/9t  =  Vm(B.  -  Br)  (B.4.3) 


Also,  there  must  be  continuity  of  the  magnetic  field  at  the  upper  boundary 
of  the  cavity,  or 

BV  =  B.+Br  (B.4.4) 
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By  combining  Equations  (B.4.3)  and(B.4.4)  to  eliminate  B^,  the  following 


differential  equation  for  Bv  is  obtained 


(H/'V  )dB  /dt  +  B  =  2B. 
m  v  vi 


which  has  the  solution  for  constant  >3  that 


(B.4.5) 


Bv  =  2Bi[>  -  “P(- 


(B.4.6) 


The  intensity  of  the  reflected  magnetic  field  at  the  boundary  at  z  =  0 
follows  from  (B.4.4)  and  (B.4.6)  as 


Br  =  B.  [l  -  2  exp(-  Vmt/H)j 


;  z  -■  0 


(B.4.7) 


The  amplitude  of  the  reflected  wave  anywhere  in  the  fluid  obviously  differs 
from  (B.4.7)  only  by  a  time  lag  factor  At  =  z/V  ,  or 


Br(z,  t)  =  B.  [  1  -  2  exp(-  (Vmt  -  z)/H)]  S(Vmt  -  z) 


for  z  >0 


(B.4.8) 


For  an  incident  wave  front  of  arbitrary  structure  where  B.(z)  is  not 
constant,  the  magnetic  field  in  the  cavity  follows  from  Duhamel’s  integral 
and  Equation  (B.4.6)  as 


B 


=  2  exp(-  t)  ^  B.  (t)  exp(+  x)  dT 


(B.4.9) 


where 


and 


x  =  V  t/H 


m 


B.(t)  =  B.(V  t/H)  for  z  =  0 
i'  r  m 
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Also  of  interest  is  the  case  of  a  steady  state  sinusoidal  descending 
wave.  For  this  case  the  magnetic  field  of  the  descending  wave  (at  z  -  0) 
can  be  represented  as 

B . (t)  -  B^  ^  exp  (iut)  ;  z  =  0  (B.4. 10) 


where 


B.  is  the  amplitude  of  the  wave, 
lm 

The  corresponding  expressions  for  the  magnetic  field  in  the  cavity  and  the 
amplitude  of  the  reflected  wave  (at  z  =  0)  are  found  to  be  of  the  form 


B  /B.  =  2/(1  +  iu>K/V  ) 
vi  m 


(B.4.11) 


B  /B.  =  (1  -  iwH/V  )/(l  +  iwH/V  ) 
r  i  rn  m 


(B.4.12) 


B.  5  REFLECTION  FROM  AN  UNLIMITED  EXPONENTIAL  BARRIER 


This  section  considers  the  problem  where  a  hydromagnetic  wave 
traveling  in  a  uniform  perfect  plasma  (with  phase  speed  VmQ)  suddenly 
encounters  a  plasma  barrier  region  of  exoonentialiy  increasing  phase  speed 


Incident  wave 


0 

Fig.  B.  3  Unlimited  Exponential  Barrier  Problem 


(at  z  =  0  in  Fig.  B,  3).  For  this  situation,  the  hydromagnetic  wave  is 
totally  reflected.  For  sinusoidal  steady  state  waves,  the  appropriate 
differential  equations  are 


whe  re 


92E/92z  +  (u2  /V^)E  =  0 
uB  =  ±  i9E/ 9z 


V 

m 


V 

mo 

V  exp(z/2p) 
mo 


;  z  <  0 

;  z  >  0 


(B.5.1) 

(B.5.2) 

(B.5.3) 
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V  is  the  initial  phase  speed 

mo 

p  is  a  scale  distance 

where  (B.5.1)  and  (B.5.2)  were  obtained  by  setting  v=  0  and  3/3t  =  iu> 
in  (B.2.19)  and(B.2.  12),  component  subscripts  being  omitted  here,  and 
where  the  upper  sign  in  (B.5.2)  refers  to  the  modified  Alfven  wave  mode  and 
the  lower  one  to  the  Alfven  mode.  The  reflection  coefficient  for  the  magnetic 
field  of  incident  waves  moving  to  the  right  in  Figure  B.3  may  be  conveniently 
obtained  with  the  aid  of  Smiley’s  (1964,11)  Equation  (3.25),  and  with  Vi  =  0  and 
written  in  a  different  form,  as 


where 


B  /B.  =  -  exp(-  2ia) 
r  i 

(B.5.5) 

a  =tan->  [j  ,(4 )  /  JQ(S )] 

(B.5.6) 

1  •  2»P/Vmo 

(B.5.7) 

B.5.2  Step- Function  Response 

It  is  next  desired  to  obtain  the  wave  form  of  the  reflected  wave  which 
is  procuced  by  the  incidence  of  a  step  function  (shock)  wave  upon  the  exponen¬ 
tial  barrier. 

It  is  first  convenient  to  replace  the  exact  expression  (B.5.5)  by  a 
more  tractable  approximate  expression.  We  choose  the  approximation 


exp(-  i|) 


£  <4  * 


B  /B.  = 
r  i 


exp  £  -  i(2£  -  j-  ir)J  ;  £  >  y-  it 


(B.5.8) 
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which  is  compared  with  the  exact  expression  (B.5.5)  in  Figure  B.4.  The 
corresponding  step-function  response  may  be  obtained  by  applying  Equations 
(E.3)  and  (E.7)  of  Appendix  E,  to  (B.5.8)  or,  more  conveniently,  by  setting 
e  -  0  in  Equation  (B.6,8)  of  Section  B.6.2  as 

-  Br  =  L  +T--1Si[^(T  -1)]  -  tt'1  Ci[|  |t  -  2]  ]  (B.5.9) 

where 

Bi  =  S(T ) 


(B.  5. 11) 


The  time  history  shape  of  the  reflected  wave,  according  to  Equation 
(B.5.9).  is  shown  in  Figure  3.18  of  Sectionlll.  It  is  seen  that  the  reflected 
wave  intensity  is  small  at  first,  then  builds  up  to  a  peak  (which  is  infinite)  in 
the  time  4p/Vmo-.  and  then  rapidly  decreases  to  a  final  reflection  coefficient 
of  minus  unity. 


B.  6  REFLECTION  AND  TRANSMISSION  AT  A  LIMITED  EXPONENTIAL  BARRIER 


B.6.1  Frequency  Response 


This  section  considers  the  problem  where  a  hydromagnetic  wave 
traveling  in  a  uniform  perfect  plasma  (with  phase  speed  VmQ)  suddenly 
encounters  a  plasma  barrier  region  where  the  phase  speed  increases 
exponentially  at  first  and  then  asympototically  approaches  a  constant  large 
phase  speed  V  (Fig.  3.  9).  This  is  the  same  problem  as  that  of  Section  B.5, 
except  that  the  phase  speed  here  is  finite  at  z  =  »,  but  was  infinite  in  Section  E.6. 

For  this  problem  the  phase  speed  is  represented  as 


Vs  = 
m 


;  z  <  0 

V'2  '  exp  {-  z/p)  +  V  "2  ;  z  >0 

m.  rni 


(B.6.1) 
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where 


is  the  initial  phase  speed 
is  the  final  phase  speed 
is  a  scale  distance 


The  differential  equations  (B.5.1)  and  (B.5.2)  remain  applicable.  The 
reflectioi  coefficient  for  the  magnetic  field  of  incident  waves  moving  to  the 
right  in  Figure  3.17  may  be  conveniently  obtained  with  the  aid  of  Smiley's 
(1964,  II)  Equation  (3.25)  as 


B  /B.  =  -  exn 
r  x 


pj+2icot-1  (vn 


«!/ J;..  (4) 


(B.6.2) 


where 


VI  * 

4  = 


(B.6.3) 


B^  is  the  amplitude  of  the  incident  magnetic  field 
B^  is  the  amplitude  of  the  reflected  magnetic  field 

We  shall  consider  only  the  special  case  of  small  V mo/Vmi  (therefore, 
small  t>j/£  ),  for  which  case  Equation  (B.6.2)  reduces  to  the  form 


-  B  /B.  = 
r  l 


exp  (-  2  i  a) 


;  I  <  N 


j^exp  j-  i(2  £  -  y  it)  -  iijvjJ]  ;  |  >: 


(B.  6. 4a) 


where  N  is  any  large  number  (>>1)  chosen  such  that  NV mQ/V =  NVj./|  is 

a  small  number  (<<  1),  and  a  is  given  by  Equation  (B.5.6).  Equation  (B.6.4) 

may  be  obtained  from  (B.6.2)  through  examination  and  application  of  the  series 

expansions  of  Erdelyi  (1953,  II)  for  Bessel  functions  with  small  and  large 

arguments,  while  noting  the  implied  restriction  here  that  v j./£  <  <  1  (since 

V  /V  <<1). 
mo  m  i 
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Another  convenient  approximation  for  Equation  (33.6. 2)  for  small 
is  the  expression 


V'5 


-Bf  /B.  =  exp  (-  2  i  a  -  it  j ) 


(B. 6.4b) 


Also  of  interest  is  the  energy  content  cf  the  transmitted  wave,  which 
may  be  obtained  from  Smiley's  (1964,  II)  Equations  (3.24)  and  (3.28)  as 


where 


St/Si  =  [4(sinh  -r  )  ]  /  J.  (|)  -  i!J.v  (t)/-J¥  4  vj  '  (B.  6.  5) 


is  the  incident  energy  flux 


St  is  the  transmitted  energy  flux 

J.  (4)  is  a  Bessel  function  of  the  first  kind  of 
1  I  imaginary  order  and  real  variable 

For  -0  and  4  large,  St/S.  may  be  evaluated  more  simply  from  Equation 

(B.6.4),  using  the  equation  S^/S^  =  1  -  J  B^/B^  j2  ,  as 


St/S^  =  1  -  exp  (-  2vvT)  for  4  large 


(B.  6.  6.  a) 


For  Vj/4  and  4  small,  another  equation  may  be  obtained  through  use  and 
simplification  oC  Smiley's  (1964,  II)  Equation  (3.71)  as 

st/si  =  4(vi^  H1  +i!2)/u  +  43/4  +  2(vI/i)(i  +  £42)];  4<i  (B.6.6b) 


Equation  (B.6.5)  is  plotted  in  Figure  3.  19  of  Section  HI  for  a  ratio  of 

V  /V  equal  to  1/30. 
mo  m  i  n 


B.6.2  Step-Function  Response 


It  is  next  desired  to  obtain  the  wave  form  of  the  reflected  wave  which 
is  produced  by  the  incidence  of  a  step  function  (shock)  wave  upon  the  exponential 
barrier.  As  in  Section  B.5,2,  it  is  convenient  here  to  first  replace  the 
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rC?  v*  ’***■ 


expression  (B.6.4)  by  a  more  tractable  approximate  expression,  which  we 
choose  as 


-B  /B.  = 
r  1 


exp  {  -  i  g)  ;  ^<7  ir 

exp  £-  i{2  £  -  £  it)  -  tr  VjJ  ;  |  >  f  rr 


(B.6.7) 


This  expression  is  analogous  to  Equation  (B.5.8)  and  is  as  accurate 

(for  Vj/  £<  <  1)  as  is  implied  by  the  comparison  of  the  exact  and  approximate 

expressions  in  Figure  B.4. 

The  step-function  reflected  wave  response  corresponding  to  (B.6.7) 
may  be  obtained  by  applying  Papoulis’  (1962)  Equation  6-21  to  the  upper  term 
in  Equation  (B.6.7)  and  by  applying  Equations  (E.  3)  and  (E.7)  of  Appendix  E 
to  tne  lower  term  in  Equation  (B.6.7).  The  resulting  equation  for  the  reflected 
wave  response  is 


Br  =  f  +  Tr_1Si[  £  { r  -  1)]  +  T  -»r{e,[^  (i  (x  -  2)  -«)]}  (B.6. 


8) 


where 

B.  =  S(t) 


t  =  Vt/2p 
00 

E,(z)  =  ^  fexp(-  z)/z  Jdz 
z 

€  =  V  /V 
mo  m  i 

x 

Si(x)  =  I  (sin  x/x)dx 


(B.6.9) 


E,(z)  is  the  exponential  integral  according  to  the  convention 
of  Abramowitz  et  ai .  (1964) 

It  may  also  be  noted  that  since  c  has  been  assumed  to  be  small  for  the 
present  problem,  the  exponential  integral  in  Equation  (B.6.9)  can  usually  be 
replaced  by  the  corresponding  expression  for  «  =  0,  which  is  (Abramowitz 
et  al . ,  1964) 
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<B. 6. 10) 


R  (Ei[t  »  -  2)]}=  -  Ciff  tr 


2 


] 


where 


co 


Ci(x)  = 


-I 


{cos  x/x)  dx 


This  will  give  a  good  approximation  except  where  Jt  -  2|  <  t  ,  near  which 
point  e  is  important. 

The  time  history  of  the  reflected  magnetic  field,  according  to 

Equation  (B.6,8),  is  essentially  the  same  as  is  given  by  the  curve  in  Figure  3.  18 

of  Section  III  for  V  =  co,  except  that  the  narrow  infinite  peak  is  reduced  to  a 

finite  value.  For  example,  for  V  /V  =  1/30,  the  peak  would  be  reduced 

mo  m  i 

from  -oo  to  the  value  -i.7. 


B.  7  REFLECTION  AND  TRANSMISSION  AT  AN  ABRUPT  BARRIER 


This  section  considers  the  reflection  and  transmission  of  hydromngnetic 
waves  .it  an  abrupt  boundary  between  two  uniform  plasma  regions  of  different 
phase  speeds  Vmi  and  Vmt  {Fig.  B.5). 


Incident  wave  i/  -  i/ 

vm  ~  vmt 


Fig.  B.5  Abrupt  Barrier  Problem 
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For  sinusoidal  steady  state  waves,  the  magnetic  fields  for  this  problem  may 
be  represented  in  the  form 


r  exp  (-  ik.  z)  +  B^  exp  (+  ik.  z)  ;  z  <0 
Et  exp  (-  ik^z)  ;  z  >0 


(B. 7. 1 ) 


where  B^,  B^,  B^  arc  the  amplitudes  of  the  incident,  reflected  and  trans¬ 
mitted  waves, respectively,  where  a  common  time  factor  exp  (itt)  has  been 
omitted  and  where 

k.  is  the  incident  wave  number 
1 

k^  is  the  transmitted  wave  number 

and  where  the  wave  numbers  are  related  to  the  corresponding  phase  speeds 
or  conductivities  by  Equations  (B.2,28)  through  (B.2.30).  Component  sub¬ 
scripts  for  the  two  possible  wave  modes  are  omitted  in  this  section  (see 
Section  B.  2). 

The  electric  field  equations  corresponding  to  Equations  (B.7.  1}  follow 
from  Equation  (B.  2.  12)  as 

E.  exp  (-  i  k.  z)  +  E  (+  i  k.  z) 

{i  r  '  1  r  1 

E^  exp  (-  i  k^  z) 

where 

E.  =  ±(<r  /  k.)  B. 

1  '  r  1 

Er  =  T(..»/  k.)  Br 
Et  =±(‘>’/kt)Br 


;  z  <0 

;  z  >0 


(B.  7.2) 


(B.7.3) 


and  where  the  upper  signs  refer  to  the  modified  Alfven  mode  and  lower  signs 
to  the  Alfven  mode. 


Application  cf  the  boundary  conditions  of  continuity  of  the  electric  and 
magnetic  fields  at  the  boundary  at  z  =  0  to  Equations  (B.7.1)  and  (B.,7.2)  gives 
the  equations 


E.  +  E  =  E 
1  r 

B.  +  B  =  B 
1  r 


*  } 


(B.7.4) 


and  combination  of  these  equations  with  Equation  (B.7.3)  gives  the  following 
reflection  and  transmission  coefficients 


Br/B.  =  -  Er/E.  =  (kt  -  k.)  /  (k,  +  k.) 

VBi  »  2  V(kt  +  ki> 

Ej/E.  =  2  k./(kt  *  k.j 


•>« 


(B.7.5) 


For  the  special  case  where  a  wave  moves  into  a  much  denser  medium 
(with  a  much  lower  phase  speed),  kfc  »k.  and  Equations  (B.7.5)  reduce  to  the 


form 


X 


B  IB.  =  -  E  / E.  =  1  -  2c, 
r  l  r  l  1 

Bt/B.  =  2  (1  -  Cj) 

E./E.  ^  2c, 
t  i  1 

ct  =  k./k  «  1 


(B.7.6) 


y 


For  the  special  case  where  a  wave  moves  into  a  much  more  rarefied 
region,  k  «  k^  and  Equation  (B.7.5)  reduce  to  the  form 


B  /B  =  -  E  /E.  =  -  (1  -  2c3)  ^  -1 
r  i  v  i  v 


Bt/B.  =  2c  2 


Et/E.  =2(1  -c; 


c3  kt/k.  «  1 


(B.7.7) 
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B. 8  ESCAPE  OF  A  MAGNETIC  PULSE  FROM  A  VACUUM  CAVITY  II 


This  section  considers  the  escape  of  a  finite  magnetic  pulse,  suddenly 
applied  to  a  vacuum  cavity,  which  is  bounded  below  by  a  perfect  conductor  and 
above  by  a  finite  uniform  plasma  region  (Region  II  in  Fig.  B.6),  which  is  in 
turn  bounded  above  by  another  region  or  set  of  regions  (Region  III  and  IV  in 
Fig.  B.6)  of  various  properties.  This  problem  is  a  generalization  of  the 
problem  in  Section  B.3  and  most  of  the  equations  derived  in  that  section  apply 
here  also  with  minor  modifications.  Also  considered  here  is  the  frequency 
response  of  the  vacuum  cavity  for  a  sinusoidal  application  of  magnetic  flux. 

B.8.1  Basic  Equations 

Equations  (B.  3.1)  through  (B.  3.  3  and  (B.  3. 6)  apply  to  the  present 
problem  for  regions  I  and  II  without  modification.  Equations  (B.3.4)  and 
(B.3. 5)  must  be  modified  to  include  the  presence  of  a  reflected  wave  moving 
downward  in  region  II,  and  may  be  generalized  as 

E  (z,  s)  =  A(s)  [exp  (-  y  z)  “  R(s)  exp  (+  y  z)]  (B.8.1) 

B'(z,s)  =  ±  A(s)  (y/s)  [exp  (-  y  z)  +  R(s)  exp  (+  y  z)]  (B.8.2) 

where 

y  =  y  (s)  =  \lss  TVs  /Vm  (B.8.3) 

and  where  R(s)  is  a  reflection  coefficient  which  may  be  defined  by  the  require¬ 
ment  that  R(ii! )  is  the  ratio  of  the  reflected  magnetic  field  to  the  incident 
magnetic  field  produced  by  a  steady-state  sinusoidal  wave  moving  upward  in 
the  plasma  region  II,  as  measured  just  above  .he  lower  edge  of  the  plasma 
region  II  (at  z  =  0+). 

The  overall  reflection  coefficient  R  may  be  expressed  in  terms  of 
reflection,  transmission  and  attenuation  factors  for  the  regions  and  boundaries 
between  regions  in  Figure  B.6  as  follows.  (It  is  convenient  here  to  think  in 
terms  of  sinusoidal  steady-state  wave  motion  with  coefficients  being  functions 
of  iut  instead  of  Laplace  transform  functions  of  s_  .)  Reflection  and  transmission 
coefficients  at  the  various  boundaries  in  Figure  B.6  are  defined  as  indicated 
by  the  various  arrow  clusters  as 


1  1  1 


Vacuum  Cavity 
(Region  I) 


Fig.  B.6  Multiple  Layer  Problem 


Rjr(ii)  is  the  reflection  coefficient  for  a  wave  crossing  the 
boundary  from  region  j  to  region  r 

1jr(i'«)  is  the  transmission  coefficient  for  a  wave  crossing 


the  boundary  from  region  j  to  region  r 


*  j(ir)  is  the  attenuation  experienced  by  a  wave  traversing 
the  region  j,  which  may  be  expressed  in  the  forms 


T]  (ii)  =  exp  [-  i  k  (ii)  h  ]  =  exp  [-  y  (ii)  h  ] 

J  J  J  J  J 


(B.8.4) 


k.(ix')  is  the  wave  number  in  region  j 


Y(i  i’)  =  i  k  (i  *  ) 


(B.8.5) 


h. 


is  the  height  of  the  j’th  region 


If  the  region  III  is  semi-infinite  (no  region  IV)  the  overall  reflection 
coefficient  R  may  be  expressed  simply  as 


2 

R  =  Rg3  r|g  (no  region  IV)  (B.8.6) 

S 

! 

j 

If  the  region  IV  is  present,  a  more  complicated  expression  results  as 
a  consequence  of  multiple  reflections  in  the  region  III.  Referring  to  the  arrow 
groups  in  Figure  B.6,  continuity  of  the  descending  magnetic  field  at  the  lower 
edge  of  region  II  (at  z  =  0+)  gives  the  relationship 


(B^  Up  Rpg 


4. 


Taj)  tv,  =  RB. 


(B.8.7) 


and  for  the  region  III,  continuity  of  the  ascending  magnetic  field  at  the  upper  edge 
of  the  region  requires 


(B.  T23  +  B 


r? 


Tfe  Ra?)  ns  =  B;. 


(B.8.8) 


where 


B  =  B .  R  *,  a 
r?  12  3 4 


(B.8.9) 
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-  < 


B^,  designate  the  magnetic  fields  of  the  ascending  and 

descending  waves,  respectively,  at  z  =  0  ’ 

Bis*  Brs  designate  the  magnetic  fields  of  the  descending  and 
ascending  waves,  respectively,  at  the  upper  edge  of 
the  region  III 

Solution  of  Equations  (B.8.6)  through  (B,8.8)  for  the  overall  reflection  factor 
R  gives  the  result 


R  ~  ’"'s  1^23  ■*"  r3  ^23  "^33  ^34  !  0  ^3  ^33  ^34  )  ] 


(B.8. 10) 


Returning  now  to  the  Laplace  transforms  of  the  differential  Equations 
(B.3.6),  (5.8.1)  and  (B.8.2),  the  unknown  functions  A(s)  and  (s)  are 
determined  by  the  requirement  of  continuity  of  E  and  B  at  z  =  0,  which  gives  the 
relationships 


A(s)[l  -R(s)]  =  TH(sB;(s)  -BvJ 


(B.8.  11) 


J\-| 

h'Via 


A  (s)  (y/s)  [l  +  R  (s)]  =  B^(s) 


(B.8. 12) 


where  upper  signs  refer  to  the  modified  Alfven  mode  and  lower  signs  to  the 

t  v 

Alfven  mode.  The  solution  of  these  equations  for  Bv(s)  is 


where 


B';  (s)  =  Bvo  (1  +  P  (S)]-1  / s 


P(s)  =  [(1  -  R)/(l  +  R>]  /  (YH) 

=  [y  H  tanh  8  ]-1 


8  =  -  &  In  (-  R) 


(B.8.  1  3) 


(B.8. 14) 


KB'- 


The  inverse  Laplace  transform  of  Equation  (B.8.  12)  gives  the  cavity 
response  to  a  step  function  flux  input.  Also  of  interest  is  the  frequency 
response  H(n)  of  the  cavity  to  a  steady-state  flux  input,  which  is  obtained  by 
multiplying  Equation  (B.8.  13)  by  s  and  replacing  s  by  ia:  ,  or 


53  f 
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H(i.)  =  1  /[  1  r  P  <i*)] 


(B.8.  15) 


where  this  expression  has  been  normalised  by  setting  B  =  1.  The  corre¬ 
sponding  characteristic  or  resonant  frequencies  occur  at  the  poles  of  H(i:), 
where 

P(i'n)  =  -1  {resonant  condition)  (B.8.16) 


is  a  characteristic  frequency 


B. 8. 2  Case  I 


The  simplest  case  of  interest  here  is  the  case  where  the  plasma 
region  above  the  vacuum  cavity  is  semi-infinite  (regions  IE  and  IV  in  Figure 
B.6  are  absent).  This  is  the  same  problem  considered  in  Section  B.  3.  For 
a  plasma  without  Hall  effects  present,  the  frequency  response  function  of  the 
cavity  for  this  case  follows  from  Equations  (B.8.  1  4)  and  (B.8.  15),  with  R  -  0, 
as 


H(i.i)  =  yH/(l  +  yH) 

=  i  k  H/  (1  +  ik  H) 


(B.8. 17) 


where  y  and  k  are  considered  functions  of  i;  .  The  corresponding  step 
function  response  is  given  by  Equations  (B.3.11),  (B.3.15)  and  (B.3.19)  for 
complex,  perfect  and  ohmic  plasma  regions,  respectively. 


Hall-effect  plasma 

For  the  case  of  a  Hall-effect  plasma,  the  step-function  response  is 
given  by  Equation  (B.  3.  34).  The  frequency  response  function  is  a  four- 
component  tensor  H„(i>  ),  whose  form  may  be  inferred  by  inspection  of 
Equations  (B.3.31)  and  (B.3.34)  to  be  as  follows 


1  .5 


H=o=(iU)'  *  A+ 


H  (iu)  =  C12  A 
xyv  1  12  - 


H  (ia)  =  Coi  A 
yx  '  21  + 


H  (id>)  =  Cm  A 
>7  ^  - 


1 

J 


A^  =  \/  ituT  ^  /  { 1  +  N/iuoT^  ) 


(B.8.18) 


(B,8.19) 


H.j(iu)}  is  the  frequency  response  in  the  i-direc&cn  due  to 
a  unit  flux  applied  in  the  j -direction 


and  where  the  C..’s  and  T..  are  defined  in  Section  B.3.5. 

B.8.  3  Case  II 


The  next  case  to  be  considered  is  that  case  where  the  region  II  in 
Figure  B.6  is  bounded  above  by  a  semi-infinite  uniform  plasma  region  III  of 
mich  lower  density  than  region  II.  In  this  case  the  reflection  coefficient  R 
is  given  by  Equation  (B.8.6)  with  Rg3  «  -1  from  Equation  (B.7.7)  of  Section  B.7, 
or 

R{s)  =  -  =  -  exp  [-2h3Y(s)]  (B.8.20) 

The  corresponding  frequency  response  function  is  (using  Eqs.  (B.8.  14), 

(B.8.  15)  and  {B.8. 16)): 


H  (iu))  =  (1  +  [yH  tanh  (yhj,)]  -1)-1 
=  (1  -  [k  H  tan  k  j"1)-! 


(B.8.21) 


and  the  corresponding  Laplace  transform  of  the  step  function  response  is 


sB^(s)  =  (.1  -r  [v  H  tanh  {yHs)!"1)-1  (B.8.22) 

In  order  to  invert  the  last  expression,  it  is  convenient  to  first  expand 
the  right  hand  s*de  cf  Equation  (B.8.22)  in  an  infinite  partial  fraction  series 
(e.g.,  Churchill,  1944),  which  leads  to  the  result 


Bv(s)  =  [l+2sCn§2/[(^-§|)s] 

(B.8.23) 

§  =  khg  =  -  iyhg 

(B.8.24) 

a  =  H/h, 

(B.8.25) 

Cn  =  a/[l+a+a»?J] 

(B.8.26) 

and  where  the  §  *  s  are  the  n  positive  roots  of  the  equation: 

a5ntan§n=1  (B.  8.27) 

B.8.3.1  Complex  Plasma 

For  the  case  of  a  complex  plasma,  Equation  (B.8.3)  applies  and 
Equation  (B.8.23)  takes  the  form 

B*(s)  =  [1  -  2  S  Cnq  /(s^+Vs  +%)]/s  (B.8,28) 

where 

?n  =  ?nVm/h=  (B.S.29) 


This  equation  may  be  inverted  with  the  aid  of  Erdelyi's  (1954,  p,  230)  trans¬ 
form  (9),  to  give  the  cavity  step  function  response  as 
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Bv(t)  =  1-22  Cn(l  “  exp  (-^yt)[cos  ui^t  +  ^(v/a)^}  sin  (B.b.30a) 


=  22  exp  (-  £  Vt)[cos  ou^t  +  £(l>/wn)  sin  u^t] 

where 


(B.  8. 30b) 

(B.8.31) 


B.8.3.2  Perfect  Plasma 

For  the  case  of  a  perfect  plasma,  Equations  (B.8.30)  and  (B.8.31)  apply 
with  v  =  0. 


B.8.3.3  Ohmic  Plasma 

For  the  case  of  a  simple  ohmic  plasma,  Equation  (B.2.30)  applies 
(s  replacing  itu),  and  Equation  (B.8.23)  takes  the  form 

-  ^ 

B;(s)  =  [1  -  22  CnTln/(s  +ry]/s  (B.8.32) 

where 

%  =  0(ki<rh2  }  -> 
and  the  inverse  transform  of  this  equation  is 


B.  (t)  =  W(a,  t/T)  (B.8.33) 

where 

T  =  |i  cr  H3  (B.8.34) 

W(a,  T )  =  1-22  C  [1  -  exp  (-  a3§3  r  )]  (B.8.35) 

=  2  2  C  exp  (-a2 1  2r) 
n  n  r  n 

and  the  £n's  are  the  positive  roots  of  Equation  (B.8.27) 
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B.8.3.4  Hall-Effect  Plasma 


Fcr  the  case  of  a  Hall-effect  plasma,  the  step  function  response 
function  may  be  directly  inferred  by  comparing  Equations  (B.3.21),  (B.3.34) 
and  (B.8.33),  with  the  result  that  the  lower  case  w’s  in  Equation  (B.3.34) 
are  replaced  by  upper  case  W’s  to  give 

Bx(t)  =  qx  W(a,  t/T+)  +  Cla  W(a,  t/TJ 
By(t)  =  Cai  W(a,  t/T+)  +  W(a,  t/TJ 


The  corresponding  frequency  response  function  is  a  four  component 
tensor  of  the  form 


where 


Hxx{iU,)  "  Cn  A+ 


H  (iu>)  =  Cla  A_ 


Hyx(i«))  =  C21  A+ 


Hyy(i«»)  =  Cgg  A_ 


(B.8.37) 


A .  =  ( 1  +  [  JicuT  tanh  ( View  T  / a)] -1 )" 1 
s:  ±  ± 


(B.8.38) 


B.8.4  Case  III 


The  next  case  to  be  considered  is  that  case  where  the  region  II  in 

Figure  B.6  is  a  uniform  perfect  plasma,  bounded  above  by  a  semi -infinite. 

region  III,  where  the  phase  speed  at  first  increases  exponentially  with 

increasing  altitude  and  then  approaches  a  constant  large  phase  speed  at  high 

altitudes.  (No  region  IV  is  assumed  present.)  The  specific  assumed  phase 

speed  variation  is  given  by  Equation  (B.6. 1),  referenced  lo  the  lower  edge 

of  Region  HI.  The  phase  speed  in  the  uppermost  part  of  region  HI 

(V  )  is  assumed  much  larger  than  that  in  the  uniform  plasma  region  II  (V  ). 
mi  °  °  mo 
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The  corresponding  overall  reflection  factor  R  for  this  problem  is 
given  by  Equation  (B.8.6)  with  given  approximately  by  Equation  (3.6.4b), 


R  =  -  exp  [-  2ia-  rr  |i/  ] 


(B.8.39) 


where 


%  =  exp  {-  ikj-hg) 


(B.8.40) 


kg  =  a1  /V  =  -  iy2 
*  mo  1  * 


(B.8.41) 


a=  tan-1  [JX(§)/J0(§)1 


(B.  5.6) 


?=  Ztop  IV  k 
vl  =  2u)p/Vt 


(B.6.5) 


is  the  phase  speed  at  the  lower  edge  of  region  III  and  in 
region  II 

is  the  phase  speed  in  the  uppermost  part  of  region  IV 

2p  is  the  scale  height  for  the  exponential  phase  speed 
variation  in  region  III. 


In  a  more  convenient  form,  Equation  (B.8.39)  may  be  expressed  as 


R  =  -  exp  (-2ia-  2ib§) 


(B.8.42) 


where 


b  =  £(hg  /pTiffe);  u;  <  0 


(B.8.43) 


<r  =  V  /V  «  1 
mo  mi 


(B. 8.44) 
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and  the  corresponding  parameter  S  {see  Eq.  (B.8.14))  is 

6=i(a  +  bs)  (B.8.46) 

The  characteristic  frequencies  for  this  system  are  given  by 
Equations  (B. 8. 14),  (B. 8.  16)  and  (B. 8. 41)  as 

i  kg  H  tanh  6+1=0 
or 

kjH  tan  (a  +  b  ?)  =  1 
or 

a  §  tan  (a+  b  §)  =  1 

where 

a  =  £  H/p  (B. 8.48) 

The  roots  of  Equation  (B.8.47)  occur  in  pairs,  each  related  by  the  equation 

5_  =  -C+  (B.8.49) 

where 

§  are  the  two  values  of  §  for  each  pair 
f  is  the  complex  conjugate  of  § 

The  frequency  response  function  for  this  system  is  given  through 
Equatiors  (B.6.3),  (B.8.14),  (B.8.15),  (B.8.42)  and  (B.8.45)  as 


(B.8.46) 

(B.8.47) 


It  is  now  convenient  to  expand  the  right  hand  sine  of  Equation  (B.8.50) 
in  an  infinite  partial  fraction  series  (e.g.,  Churchill,  1944),  which  gives  the 
result 

H(iui)  =  1  +  2  A  /(?  -  §  ) 
n  n  n' 

where 

A  =  a?  /[a  +  g  (1  +  a2?s)l 
n  n  1  6nx  ’a'- 

gn  =  b  +  da/d?  for  ?  =  ^ 

=  b  +  i  -jo(§n)j,(?n)/[f 

and  where  the  §  ’  s  are  the  roots  of  Equation  (B.8.47).  This  series  may  be 
simplified  by  considering  the  symmetry  of  this  problem  (e.g.,  see  Eq.  (B.8.49)) 
as 


(B.8.51) 


(B.8.52) 


>n  ^n* 


(B.8.53) 


H(io  )  =  1  +  S  An/(§  -  §n)  -  2  An/(?  +  ?n) 


(B.8.54) 


and  where  the  summation  here  is  only  for  the  (nearly)  real  and  positive  roots 
of  Equation  (B.8.47). 

Finally,  the  step  function  response  may  be  conveniently  obtained  by 

£*€ 

replacing  H(itu)  in  Equation  (B.8.54)  by  sB  (s),  and  im  by  -  is,  and  by 

Sje 

inverting  the  resulting  function  B  (s),  considered  as  a  Laplace  transform. 
The  resulting  step  function  response  is 


Bv(t)  =  1  -  2R(2  Cm  [1  -  exp  (iu)Qt)]) 


(B. 8. 55a) 


=  2R2C  exp  (iu)  t) 
m  r  '  n  ' 


(B. 8. 55b) 
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r;r  P5C-'. 


f  a 


APPENDIX  c 


VARIOUS  INTEGRALS 


The  following  integrals,  not  listed  in  elementary  tables  of  integrals, 
were  of  interest  for  the  present  studies 


CaJx  cos  x  dx  =  +  \fx  cin  x  -  n/it/2  S°  (\/ 2'x.j-) 
C'v/jT sin  x  dx  =  -  s/x~ cos  x  +  -rr/2  C  (\/ 2x/ir  ) 


(C.  1) 


(C.2) 


where  S  (x)  and  C  (x)  are  the  Fresnel  integrals: 


A 

C  (x)  =  ^*  cos  (rx2/2)  dx 


S  (x)  =  C  sin  (irx2/2)  dx 


(C.3) 


(C.4) 


The  following  two  integrals  were  obtained  by  Foster  (Wait  and  Ccmpbell,  19**'' 


pfg.p]  =  1  exP  ("  u  g)  JQ  ftp)  u  X  d  X 


=  exp  (  -  y  s'g2  +  p2  )N gz  +  p2 
00 

J[S,p]  =  ^  exP  ("  u  g)  Jo  (XP)  u“‘  d  x 
o 

=  v(Vg2  +  P2  -  g)J  Y(V; 


(C.5) 


(C.  6) 


(C.  7) 


where 


u  =  'JxTTvr~ 


g2  +  P2  +  g)i  (C.  8) 


(C.  9) 
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From  Erdelyi  (1954,  G,  8.4  (10)) 


Q [g»  p] 


00 

=  \  exp  j-ug)  p)  «i~l  dX  (C.  10) 

J 
o 


u 


(C.  11) 
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APPENDIX  D 


EVALUATION  OF  INTEGRALS 


This  appendix  c-msiders  the  evaluation  of  integrals  of  the  form 


where 


y  =  J £{*)  gnfr) 


g  {x)  =  sin  x  -* 

1  } 
g  (x)  =  COS  X  J 


(D.l) 


(D.2) 


where  f  (x)  is  primarily  a  slowly  varying  function  inthe  sense  that  f  (x)  is 

relatively  constant  between  adjacent  zeros  of  g^{x).  hi  addition,  effects  of 

step  changes  in  f(x)  and  variations  of  the  form  f(x)  =  *>/x  a.re  considered. 

Considering  g  (x)  first,  evaluation  of  Equation  (D.  1)  over  any  range 
«  1 

of  x  where  f  {x)  is  constant  gives 


y  =  J  |f(a)  r  f  (x)(x  -  a)  J  sin  x  dx 
a 

and  using  Pierce  and  Foster's  (1956)  Equation  345 

y  =  ( j^-i(a)  +  a  f  (x)J  cos  x  +  f  (x)  sin  x  -  f  (x)  x  cos  x) 


(D.3) 


,-b 


■a. 


or,  in  a  more  convenient  form 

y  =  ["-  f(x)  cos  x  +  f  (x)  sin  x^ 

l  L  -> 


;  f  (x)  =  constant 


(D.4) 


It  is  now  assumed  that  i{x)  can  be  represented  adequately  in  the  range 
a  <  x  <  b  by  a  series  of  straight  line  segments  with  slope  changes  at  the 
zeros  of  sin  x,  as  indicated  by  the  solid  line  in  the  sketch  below.  (That  is. 
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it  is  assumed  that  the  integral  j  f{x)  dx  can  he  evaluated  adequately 


•  's  Denote  zeros  of  g(x) 


by  the  trapezoidal  rule  with  intervals  equal  to  the  half- period  of 
gn(x).)  It  then  follows  through  progressive  application  of  (D.4)  to  each 
segment  in  the  sketch  that 


=  :(x)  cos  x  T  f  (x)  sin  xj 


which  is  the  same  form  as  (D.4),  but  without  the  restriction  that  f’(x)  be 
constant  over  the  entire  interval  a<  x  <  b. 

The  corresponding  expressions  for  y2  maybe  obtained  similarly 
through  use  of  Pierce  and  Foster's  (1956)  Equation  349  as 


y2=|f(x)  sin  x  +  f*(x)  cos  xj 


If  the  function  f(x)  contains  discontinuities  in  slope  which  do  not 
occur  at  the  zeros  of  sin  x  or  cos  x,  it  is  evident  from  inspection  of 
Equations  (D.4)  through  (D. 6)  that  there  are  then  additional  contributions 
Ayj  and  Ay2  to  the  integrals  Yi  and  yz  of  the  form 


Ayj  =  -  2  6  f  (x  )  sin  x 
n  1  n  n 


Ay2  =  -  2  5  f  (x  )  cos  x 
n  n  n 
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where 


xn  is  the  value  of  x  at  the  n’th  discontinuity  in  f{x) 

,*  t 

61  (xn)  is  the  increase  (jump)  in  f  (x^)  at  x^  with  increasing  x 

An  additional  generalization  may  be  obtained  as  follows  for  the 

case  wnere  f(x)  behaves  as  required  above  except  near  a  point  x  where 

m 

it  behaves  as 


fj(x)  =  A  \/x  -  X  ;  x  —  X 

m  m 


(D.9) 


where  A  is  a  slowly  varying  function  of  x  and  f*,(x)  is  seen  to  be  infinite  at 

x  =  xm.  The  effect  of  this  infinity  of  f  (x)  on  the  integral  (D.  5)  may  be 

obtained  by  integrating  (D.  1)  between  the  limits  x,=x  -  6,  and  x,  =  x  +6, 

m  *  m 

where  6j  and  62  are  large  constants  and  A  is  considered  constant.  With  the 
aid  of  Equations  (G.  1)  and  (C.2)  of  Appendix  Q  this  integral  is  obtained  as 


sin  x  dx 


f  j(x)  COS  X 


+  7^72  A  jcos  xm  C  (72 (x  -  xm)/ir)  -  sin  xnS"(-JZ(x  -  xm) /tt)  1 1  (£>.  10) 

J  xi 

ihe  Fresnel  integrals  C  and  S  may  be  evaluated  from  the  first  two  terms 
in  their  asymptotic  series  expansion  for  large  real  arguments  (Martz,  1964, 

Eq.  (11))  as 


a. 

C  (x)  =  &  +  sin  (£tt  x?)/(irx) 

a. 

s  (x)  =£-  cos  (£  IT  x2)/(trx) 


X  >>1  and  real 


(D.  11) 


and  for  imaginary  arguments  by  the  relationships  (Abramowitz,  et.  al. ,  1964) 


C^(ix)  =  iC^(x)  j 

S^(ix)  =  -  iS"(x)  J 


(D.  12) 
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It  is  assumed  that  the  function  f(x)  is  everywhere  slowly  varying,  except 
near  x  ,  where  it  va.ries  according  to  Equation  (D.  9),  and  that  f(x)  vanishes 
exponentially  at  infinity.  It  then  follows  from.  Equations  (D.  5),  (D.  6),  (D.  14), 
(D.15),  (D.  17)  and  {D.  18)  that  the  integrals  (D.  I6)nave  the  values 
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f(x)  JQ  (x)  dx  =  -  f(x1)J1r(x1) 


Vxj  fd/dx  [  f(x)  A/x  ]) 


Jo 

{X  =  Xj. 


oo 

f{x)  Ji  (x)  dx  =  £(xi)J^(xi) 

Vxx  (d/dx  [f(x)/7x  ])  J  +‘fn2  A"  exp  (-  ix  ) 


APPENDIX  E 


FREQUENCY- TIME  TRANSFORMS 


This  appendix  is  concerned  with  the  development  and  tabulation  of  the 
impulse -response  and  step-function-response  of, linear  systems  characterized 
by  various  frequency  response  functions.  The  following  conventional  defini¬ 
tions  are  used. 

If  a  steady-state  sinusoidal  input  of  the  form  exp(itot)  is  applied  to  a 
linear  system,  the  output  will  be  of  the  form  H(iio)exp(L>t)  where  H(ito)  is  called 
the  complex  frequency  response  or  system  function  or  transfer  function  of  the 
linear  system. 

The  corresponding  response  h(t)  of  a  linear  system  to  a  unit  impulsive 
input  6(t)  is  given  by  the  following  equivalent  Fourier  transform  equations 
(Papoulis,  1962) 


cc 


where 


h(t)  =  (2-rr )  1  J  exp  (itot)  dto 


-oo 


oo 


- 1 


/[’ 

5  l- 


=  tt  /  J  R(gj)  cos  tot  -  X(to)  sin  cot  j  dco 
0 


'] 


“  00 

=  ir_1  Re  /  H(i 
-  o 


H(ito)  exp  (itot)  dto 


(E.l) 


(E.2) 


(E.3) 


R(co)  is  the  real  part  of  H(ito) 

X(co)  is  the  imaginary  part  of  H(ito) 

Re£x  J  designates  the  real  part  of  £x  J 
h{t)  is  the  response  to  a  unit  impulse  6 (t)  where  h(t)  is 
assumed  to  be  a  real  function 
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and  for  the  simpler  case  of  causal  linear  systems  (Papoulis,  1962),  where 

•A# 

the  response  h(t)  is  zero  for  negative  time,  the  following  simpler  equations 
apply 


oc 

h(t)  =  (2/ir)  J  R{w)  cos  wt  dw 

(E.4) 

0 

CO 

=  -  (2/tr)  J  X(w)  sin  wt  dw 

(E.5) 

0 

=  L-1  [h(s)] 

(E.6) 

where  L-1  [f (s)]  is  the  inverse  Laplace  transform  of  f(s),  subject  to  the 
restriction  that  f(s)  is  analytical  in  the  right  half  of  the  complex  s -plane 
where  the  real  part  of  s  is  positive. 

The  response  of  a  linear  system  A(t)  to  a  unit  step  function  S(t) 
is  in  general  given  by 


where 


t 

A(t)  =  J  h(t)  dt  (E.  7) 

-oo 


A(t)  is  the  response  to  a  unit  step  function  S(t) 


and,  for  the  case  of  causal  systems,  by(Papoulis,  1962) 


oo 

'*)  J  £r(u>)/wJ 


A(t)  =  (2/ir)  (  R(w)/w  sin  wt  dw 


=  R{0)  +  (2/  ir  )  J  £  X(w)/w  j 
=  L~ 1  [H(s)/s] 


cos  wt  dto 


(E.8) 

(E.9) 

(E.  10) 


While  real  physical  systems  are,  of  course,  all  causal  (i.e.,  output  does 
not  precede  input),  the  approximations  involved  in  deriving  some  equations 
in  this  report  result  in  a  few  non -causal  expressions. 
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Table  E.l  presents  expressions  for  h(t)  and  A(t)  computed  according 
to  the  above  equations  for  a  variety  of  frequency  response  functions  of  interest 
in  other  parts  of  this  report.  The  remainder  of  this  section  consists  of  a 
partial  listing  of  the  integral  transforms  used  to  determine  each  listed  function. 

Item  1.1:  Inverse  Laplace  Transform;  Churchill  (1944),  Appendix  II  (14); 

Erdelyietal.  (1954,  I),  5.6  (35);  Papoulis  (1962),  Eq.  (3-4). 

Item  1.2;  Papoulis  (1962),  Eq.  (3-4). 

Item  1.3:  Inverse  Laplace  transform;  Churchill  (1944),  Appendix  III  (82,83). 

Item  2,1a:  Inverse  Laplace  transform;  Churchill  (1944),  Appendix  11(14); 

Erdelyi,  et.  al.  (1954,  I),  5.6  (36)  and  (37). 

Item 2.1b;  Differentiation  of  2.1a  with  respect  to  c. 

Item  3.2:  Inverse  Laplace  transform;  Churchill,  Appendix  III  (8,  12). 

Item  3.3:  Inverse  Laplace  transform;  Churchill,  Appendices  II  (14)and 
III  (37,43). 

Item  5:  Papoulis  (1962),  Eq.  (3-3). 

Items  6  and  7:  Papoulis  (1962),  Eq.(3-32). 

Item  8.1:  h(t)  determined  by  convolution  of  Items  1  and  4. 

Item  8. 2:  Inverse  Laplace  transform;  Churchill  (1944),  Appendices  II  (12,  13) 
and  TJX  (8). 

Item  8.3:  Inverse  Laplace  transform;  Churchill  (1944),  Appendices  II  (14) 
and  IH  (87). 

Item  9:  Equations  (E.3)and  (E.7);  Grobner  and  Hofreiter  (1958),  p.  128  (67a). 

Item  11:  Equations  (E.3)  and  (E.7);  Grobner  and  Hofreiter  (1958),  p.  127  (66b). 

Item  13:  Equations  (E.3)  and  (E.7);  Erdelyi  et  al.  (1954,  I).  1.9  (1);  Dwight 

(1947),  679.10. 

Item  15:  Take  imaginary  part  of  expression;  apply  Equation  (E  .  5),  using 
Erdelyi  et  al  (1954,  I),  2.12(27). 
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*Sce  also  2.  la  and  2.  lb  on  next  page 


Table  E  .  1  (Cont) 
Table  of  Transforms 


Item 

H  (iw) 

A  (t) 

2.1a 

V  exp  (-  cy) 

c-1f  exp  (-  x)  J^6{x  -  f)  +  lQtiJ  -  f2)S  (x  -  f)  +  xli  (t/x2  -  f2)  j 

2.1b 

y2  exp  (-  cy) 

c"  2f 2exp  (-  x)  j^6i  (x  -  f)  +  (1  +  5  f)  6  (x  -  f) 

+  (f  x/fx2  -  f2))Io(Aa  -  f2)S(x  -  f)  +  f(l  -  2  x/ (x2  -  f2))4  (Vx2-!2)] 

I?(x,f)  =  L  (/x2  -  fa)  S  (x  -  f)A/xa  -  fs 


(E.  11) 


U(x,  f)  -  exp  (-f)  S  (x  -  f)  +  f 


/  * 


(x,  f)  exp  (-x)  dx  (E.  12) 


The  step  function  S (t) ,  the  delta  or  impulse  function  6(t)  or  SQ(t), 
and  the  n’th  derivative  of  the  delta  function  8^(t)  are  defined  as  (Papoulis, 
1962) 


f  .)  ;  t  <  0 

S(t)  =  < 

l  l  ;  t  >  o 

(E.  13) 

8(t)=  6Q(t)  =  dS(t)/dt 

(E.  14) 

6  (t)  =  dn5/dtn 
n 

(E.  15) 

oo 

j  6n(t  -  to)  f(t)  dt  =  (-  l)n(dn/dtn)  f(to) 

(E.  16) 

•  oo 
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